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ABSTRACT 

We explore explosions of massive stars, which are triggered via the quark-hadron phase transition 
during the early post bounce phase of core-collapse supernovae. We construct a quark equation of 
state, based on the bag model for strange quark matter. The transition between the hadronic and 
the quark phases is constructed applying Gibbs conditions. The resulting quark-hadron hybrid equa- 
tions of state are used in core-collapse supernova simulations, based on general relativistic radiation 
hydrodynamics and three flavor Boltzmann neutrino transport in spherical symmetry. The formation 
of a mixed phase reduces the adiabatic index, which induces the gravitational collapse of the central 
protoncutron star. The collapse halts in the pure quark phase, where the adiabatic index increases. 
A strong accretion shock forms, which propagates towards the protoncutron star surface. Due to the 
density decrease of several orders of magnitude, the accretion shock turns into a dynamic shock with 
matter outflow. This moment defines the onset of the explosion in supernova models that allow for 
a quark-hadron phase transition, where otherwise no explosions could be obtained. The shock prop- 
agation across the neutrinospheres releases a burst of neutrinos. This serves as a strong observable 
identification for the structural reconfiguration of the stellar core. The ejected matter expands on a 
short timescale and remains neutron-rich. These conditions might be suitable for the production of 
heavy elements via the r-process. The neutron-rich material is followed by proton-rich neutrino-driven 
ejecta in the later cooling phase of the protoncutron star where the z^p-process might occur. 
Subject headings: equation of state - dense matter - gravitation - shock waves - supernovae ~ neutron 
stars 



1. INTRODUCTION 

The explosion mechanism of massive stars is an active 
subject of research in theoretical astrophysics. In this 
article, we explore a scenario where the explosion is trig- 
gered due to the phase transition from hadronic matter 
to quark matter. We present specific and distinguished 
characteristics for possible explosion observables. 

Stars with an initial main sequence mass above 10 Mq 
produce extended Fe-cores, as the end product of nuclear 
burning. During the following evolution, the photodis- 
integration of heavy nuclei as well as electron captures 
reduce the pressure of the Fe-core. The initial contrac- 
tion of the core proceeds into a collapse, during which 
density and temperature rise. At nuclear densities the 
repulsive nuclear interaction causes a significant stiffen- 
ing of the equation of state (EoS), which in turn leads to 
the formation of a shock wave. As it propagates outward 
it continuously looses energy due to the dissociation of 
in-falling heavy nuclei from the progenitor. The central 
object formed at core bounce is a hot and Icpton-rich 
protoneutron star (PNS). The shock propagation across 
the neutrinospheres, i.e. the neutrino energy and fiavor 
dependent spheres of last scattering, releases a burst of 
electron neutrinos - known as the deleptonization burst - 
emitted from electron capture at free protons. The elec- 
tron neutrino luminosity rises up to several 10^^ erg/s 
(depending on the progenitor) on a timescale of 5-20 ms 



after bounce. This enormous energy loss, in combina- 
tion with the dissociation of the in-falling heavy nu- 
clei, turns the dynamic shock into a standing accretion 
sho ck (SAS) already at about 5 ms after bounce (see 
e.g. IHillebrandt fc Miillei< [l98T|: IMavle fc WilsonI [19871: 



lBruennlll989l : lMvra fc BludmanHl989f ). The deleptoniza- 
tion near the neutrinospheres results in a low proton-to- 
baryon ratio (given by the electron fraction) of = 0.1- 
0.2. The PNS interior, which did not experience shock 
heating, stays slightly less neutron-rich where Ye — 0.3 
at bounce (depending on the equation of state and the 
electron capture rates). 

In an attempt to explain core-collapse supernova explo- 
sions, neutrino heating wa s suggested as a mechan ism for 
the revival of the SAS bv lBethe fc WilsonI ([T985I) . lead- 
ing to so called neutrino-driven explosions. In spherical 
symmetry, this mechanism has been shown to produce 
explosions only for the low mass 8 .8 O-Ne-Mg-core 
fro m Nomoto (1983,1984 1987) by IKitaura et all (|2006l ) 
and IFischer et al.l ()2010bD . Multi-dimensional phenom- 
ena, such as rotation, convection and the development 
of fluid instabilities, have been sho wn to increase the 
neutrino he ating e fficiency (see e.g. iMiller et al.l [l993l: 
Herant et al . 1994: Burrows et al.lll995HJanka fc Mueller] 
1996( ). Such models help to aid t he understanding of 
aspherical explosion s (see also iBruenn et al.l 120091 : 
iMarek fc Jankal [20091 and references therein). Besides 



2 



Fischer ct al. 



neutrino heating (jBethe fc WilsonI |1985[ ). alternative 
explosive scenarios have been expl ored, such as the 
dumping of magnetic energy by iLeBlanc fc WilsonI 
(I1970D. iMoiseenko fc Bisnovatvi-Koganl ()2007f i and 
Takiwaki ct al.' f2009f) as well as by acoustic energy by 



Burrows ct al- (,20061 ). The later two scenarios require 



multi-dimensional models. 

An additional scenario is obtained by looking at the 
conditions achieved at the PNS interior during the post 
bounce phase of core-collapse supernovae. Due to the 
mass accretion on the order of several 0.1 Mq/s, the 
central density rises above several times nuclear matter 
density and the temperature above several tens of McV 
(depending on the progenitor model). With increasing 
degeneracy, the central electron fraction reduces below 
Ye < 0.25 on timescales on the order of several 100 ms 
post bounce. The timescale for the PNS contraction de- 
pends on the mass accretion rate which in turn is in 
direct correlation with the progenitor model. The con- 
ditions obtained in the PNS interiors raise the question 
about the state of matter and possibly the appearance 
of exotic matter such as hyperons and quarks. 

Assuming a purely hadronic composition, matter is 
composed of color charged quarks which cannot prop- 
agate freely and are confined into color neutral baryons 
and mesons. At high temperatures or baryon densities, 
hadronic matter is expected to change its state to a phase 
of deconfined quarks and gluons, where chiral symmetry 
is restored. Thereby, deconfinement and asymptotic free- 
dom allow color charged quarks and gluons to move as 
free and almost non-interacting particles. Due to chiral 
symmetry restoration, the quarks obtain their current 
mass values, which for the up and down quarks are in 
the range of few MeV, while the strange quark current 
mass is about 100 MeV. 

In Quantum Chromodynamics (QCD), which is be- 
lieved to be the fundamental theory behind the strong 
force, the interaction of quarks is characterized by an 
effective coupling constant a^, which monotonically ap- 
proaches zero at large momentum transfer. For small 
values of a^, a perturbativc treatment of the strong in- 
teraction can be applied. However, at the energy scale of 
quark deconfinement, the interaction strength becomes 
too large and a perturbativc treatment is no longer ap- 
plicable, requiring alternative approaches. One possi- 
ble descr iption is th e phcnomcnological quark bag model 
(see e.g. iDetar &:"T)o noghuc 198j), based on the spatial 
confinement of quarks. In Nambu-Jona-Lasinio models, 
quarks are implemented as quantum fields with point- 
like interactio ns in a chirally symmetric Lagrangian (see 
iBuballal I2005L and references therein). More advanced 
methods to describe quark matte r are based o n e.g. 
the Dyson-Schwingcr equation (see iFisched 120061 ) or ab 
initio calculations of QCD on the lattice. Ab initio 
calculations show a smooth crossover from a hadronic 
phase to a deconfined quark gluon plasma phase, at 
hig h temperatures and low baryochemical potentials (see 
e.g. iFodor fc Katzll2004l : lAoki et al.ll2006l ). On the other 
hand, field theoretical and phenomenological calculations 
point to a first order phase transition from hadronic 
to quark matter at large chemical potentia l s and finite 
temperatures (seelBerges fc Raiagopallll999t iKlevanskvl 
[T99a IPisarski fc Wilczeklll984D The many open ques- 
tions concerning critical temperatures and densities, the 



order of the phase transition and whether both chiral 
symmetry restoration and deconfinement happen simul- 
taneously or not, make QCD phase transitions an active 
topic of experimental research at the Rclativistic Heavy 
Ion Collider (RHIC) in Brookhaven, the Large Hadron 
Collider (LHC) at CERN and the Facility for Antipro- 
ton and Ion Research (FAIR) at GSI, Darmstadt and 
at NICA in Dubna (Russia). The experiments at the 
LHC and RHIC are focused on the production of matter 
at high temperatures and low baryochemical potentials. 
On the other hand, FAIR and NICA will be designed to 
study nuclear matter at large densities and temperatures, 
which are conditions that correspond to core-collapse su- 
pernova interiors. They will enable a crosscheck of results 
from terrestrial experiments and astrophysical observa- 
tions. However, the conditions in supernovae, neutron- 
rich matter and the possible production of strange quarks 
by weak interactions, differ from those in heavy-ion col- 
lision experiments. As will be discussed below, a quark- 
hadron phase transition at large baryon number densi- 
ties in heavy-ion collisions can be compatible with low 
critical densities for the appearance of quark matter in 
supernovae and neutron stars. A detailed discussion of 
this issue will be given in § 2. 

The first attempt to investigate the quark-hadron 
phase tr ansition in simul a tions of massive stars goes 
back to iTakahara fc Sato! (I1988D . They tried to find 
a connection between the appearance of quark matter 
and the apparent but statistically insigni ficant multi- 
peaked neutrino signal f rom SN1987a (see iHirata et al.l 
Il988t iBionta et al.lll987l ). using general rclativistic hy- 
drodynamics in spherical symmetry and a polytropic 
EoS. More sophistic ated input physics was applied in 
the investigation bv iGentile et al.l ()1993[ ). Their model 
was based on general relativistic hydrodynamics and a 
parametrized description of the deleptonization during 
the Fc-corc collapse as well as an improved EoS for both 
hadronic and quark matter. The transition between 
these two phases was constructed via a co-existence re- 
gion, i.e. a quark-hadron mixed phase. They found a 
significant softening of the EoS in the mixed phase as 
well as the formation of a second shock wave as a direct 
consequence of the quark-hadron phase transition. The 
secondary shock followed and finally merged with the 
Fe-core bounce shock, only milliseconds after its appear- 
ance. However, both studies were not able to predict any 
possible features in the neutrino signal emitted due to 
the lack of neutrino transport. Simulations of the QCD 
phase transition of very massiv e (100 M(^) progenitors 
hav e been reported recentl y by iNakazato et al.l ()2008aD 
and lNakazato et all (l2010al ). They applied the quark bag 
model for the quark matter EoS with a large bag con- 
stant and hence a high critical density for the onset of 
deconfinement. The transition from hadronic matter to 
quark matter was modeled via an extended mixed phase 
applying the Gibbs conditions. They confirmed a signif- 
icant softening of the EoS in the mixed phase due to the 
reduced adiabatic index. A characteristic neutrino sig- 
nature was not found due to the immediate formation of 
a black hole during the phase transition. However, they 
observed a significant shortening of the post bounce ac- 
cretion time until black hole formation, due to the soft- 
ening of the EoS at conditions where quark matter was 
found to occur. 
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TABLE 1 

List of neutrino reactions considered, where v - 



{Ue,l^e,iy^/r,'^^i/T} AND N = {n,p}. 



Label 



Neutrino reaction 



Reference 



- p ^ n + Ue 

- n ^ p + Ue 

(A,Z) ^ (A,Z-l) + Ue 

u + N ^u' + N 
V{A,Z)^u' + {A,Z) 
± ^ ,,' j_ c± 



V 



Bruenn (19851 
Bruenn (1985) 
Bruenn (1985) 
^ruenn (1985) , ^ezzacagga^^^ruenn (1993b) 
^ruenn (1985) , ^ezza£a£g^_&_^ruenn (1993a) 
^ruenn (1985) , ^ezza£a£g^^_^ruenn (1993c) 
Bruenn (1985 ). Mezzacappa fc Messer (1999') 
Thompson fc Burrows (2001') 
Buras et aL (2003), Fisch er et aL (2009') 



The present artic le foUows a d i fferent approach. It is 
the continuation of ISagert et al.l (j2009() . where first re- 
sults were discussed. We construct a quark- hadron hy- 
brid EoS, i.e. using the microscopic EoS from lShen et al.l 
(fT998i) for hadr ons and the bag model based EoS for 
strange quark matter. The resulting EoS is applied to 
core-collapse supernova simulations of low and interme- 
diate mass Fe-core progenitors. Our model is based on 
general relativistic radiation hydrodynamics and three 
flavor Boltzmann neutrino transport in spherical sym- 
metry. For these progenitors, spherically symmetric sim- 
ulations do not result in explosio ns via the delayed neu- 
trino heating mechanism (see e.g. Rampp fc Janka 1120001: 



Liebendorfer et al.ll2001bl : iJanka ll200iriThompson et al.l 
200.-^ . Here, we explore the possibility that a quark- 



hadron phase transition during the early post bounce 
evolution can result in an explosion and discuss the ob- 
servational consequences. The temperatures on the order 
of tens of MeV and the large neutron excess of matter 
in supernova interiors, enables a phase transition at low 
critical densities close to nuclear saturation density (de- 
pending on the bag model parameters). In this respect, 
the physical conditions in core-collapse supernova interi- 
ors differ from matter in heavy-ion collision experiments, 
where temperatures on the order of few 100 MeV are ob- 
tained but matter stays more or less isospin symmetric 
with an electron fraction of Ye — 0.5. The expected crit- 
ical density is on the order of several times nuclear sat- 
uration density in heavy-ion collision experiments. Fur- 
thermore, the supernova timescales (1-100 ms) and the 
timescale in heavy-ion collisions (~ 10~^^ seconds) are 
substantially different. These differences will be further 
discussed in §2. In the core-collapse supernova simu- 
lations under investigation, the phase transition takes 
place within the first 500 ms post bounce (depending on 
the critical conditions and the progenitor model). We 
observe an accelerated PNS contraction due to the soft- 
ening of the EoS in the mixed phase, in comparison to 
simulations using purely hadronic EoSs. The contraction 
proceeds into an adiabatic collapse when the maximum 
stable mass, depending on the hybrid EoS, is reached. 
During the collapse, density and temperature increase 
and a pure quark phase develops where the EoS is signif- 
icantly stiffer than the EoS in the mixed phase. The col- 
lapse halts and a strong second shock front forms. This 
shock wave accelerates at the PNS surface where the den- 
sity drops over several orders of magnitude, resulting in 
explosions with matter at outflow velocities on the order 
of 10^ km/s. Furthermore, the shock propagation across 
the neutrinospheres which are located in the hadronic 



phase, releases an additional neutrino burst dominated 
by j^e and (f^/r, i^/i/r)- We analyze the emitted neutrino 
signal from this explosion scenario for several examples 
of massive progenitors. The possible detection of the 
QCD induced neutrin o burst has been investigated by 
iDasgupta et al.l (I2010D f or a G alactic event, based on re- 
sults from ISagert et al.l (|2009f ). 

The manuscript is organized as follows. In § 2 we 
briefly discuss our core-collapse supernova model includ- 
ing the standard hadronic EoS. Furthermore, we explain 
the construction of the hybrid EoS in detail. We ex- 
plore different choices of parameters and discuss the cor- 
responding phase diagrams and hybrid star mass-radius 
relations. § 3 is devoted to the illustratio n of the results 
obtained for the 10.8 progenitor from lWooslev et aD 
(|2002f l. that will be used as the reference model. The 
results are discussed in § 4, where we compare the ref- 
erence model with simulations performed with different 
progenitor masses and different hybrid EoSs. Further- 
more, we explore the nucleosynthesis relevant conditions 
and draw conclusions about the possibility of reaching 
high magnetic fields on the order of 10^'^ — 10^^ G dur- 
ing the PNS evolution with quark matter. We close the 
manuscript with the summary in § 5. 

2. NEUTRINO RADIATION HYDRODYNAMICS 

The spherically symmetric radiation hydrodynamics 
model AGILE-BOLTZTRAN is based on three flavor 
Boltzmann neutrino transport. It solves the evolution of 
the neutrino distribution functions j^, (phase-space func- 
tions) for each neutrino flavor v = {z/g, Pe, ^jilr^ ^ti/T\^ 



dt 



dt 



(1) 



collisi( 



depending on the phase-space derivative (l.h.s of Eq.([T])) 
and neutrino- matter interactions (r.h.s of Eq.(IT])). The 
neutrino reactions for the various interactions considered 
are listed in Table [1] including the references. In spher- 
ical symmetry, the distribution functions depend on the 
phase-space coordinates system time t, baryon mass a, 
the cosine of the neutrino propagation angle fi = cos 9 
and the neutrino energy E, f,^{t,a, fi, E), where the lat- 
ter two are the momentum space representations. The 
original Newtonian va riant of AGILE-BOLTZTRAN (see 
Mezzacappa fc Bruen n 1993330) bas been extended by 
Liebendorfer et aLT i 2001al ) to solve the general relativis- 
tic equations, based on non-stationary spacetime repre- 
sented by the following line element 



ds' 



-a{t,afdf + 



r{t,a)' 
r(t,a) 
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Baryon density, log,Q(P [g/cm ]) 
8 9 10 11 12 13 14 15 
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Baryon density, log,Q(p [g/cm ]) 
6 7 8 9 10 11 12 13 14 15 
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Fig. 1. — Phase space covered in core-collap se super nova simulations at the example of two progenitor models from IWooslev &: Weaver! 
II1995I ). both using the pure hadronic EoS from lShen ct al. (1998). The white background corresponds to densities and temperatures which 
are not obtained, where the colored domain belongs to densities and temperatures which are obtained. Color-scaled is the electron fraction 
which ranges from neutron matter (blue) Ye = to symmetric matter (red) Ye ~ 0.5. 



where dfl = dO^ + siv? Odcjx^ is the hne element of a 2- 
sphere of radius r in spherical coordinates {9,(f>), and 
r' = dr/da. Note that the metric functions a and F de- 
pend on the coordinates system time t and baryon mass 
a. The evolution equations for energy and momentum 
arc obtained from the covariant derivative of the stress- 
energy tensor, ViT''^ = 0, with the following choice of 
the stress-energy tensor 

T" = p(l + e + J), T'"'=p + pK, T"* ^T*" ^ pH, 

T'^^^T'^'^ ^p+^p{J -K), 

where p, p and e are the rest mass density, matter pres- 
sure and the internal energy density. J, H and K are 
the neutrino moments, which are given by the momen- 
tum integr als of the neutrino distrib ution functions (for 
details, see iLiebendorfer et al.11200'3 ). Special attention 
has been devoted to accurately conserve lepton num- 
ber, energy and nr o ment um, as described in details in 
ILiebendorfer et all (|2004[ ). General relativistic effects 
have been investigated with respect to the PNS collapse 
and black hole formation for several massive progenitor 
stars of 40 and 50 Mq in lFischer et alj (|2OO90 . 

The numerical implementation of the above men- 
tioned physical equations, including the finite differenc- 
ing mass conserving sch eme . has been introduced i n 
ILiebendorfer et all (|2002| ) and lLiebendorfer et all (|2004f ). 
where recent improvements of the finite differencing 
scheme for the evolution of the enclos e d bary on mass 
have been introduced in iFischer et al.l ()2010b[ ). With 
these improvements, stable solutions for the evolution 
equations can be obtained for long simulation times on 
the order of seconds. 

In the following subsections, we will introd uce the EoS 
where the sophisticated baryon EoS from iShen et ahl 
(|1998D has been extended with a quark-hadron mixed 
and a pure quark phase at high densities and tempera- 
tures. 



2.1. The EoS in core-collapse supernova simulations 

The EoS has to handle a large variety of physical 
conditions obtained in core-collapse supernova simula- 
tions. The covered domain is illustrated in Fig. [1] for 
the exa mples of the 15 an d 40 M p, progenitor mod- 
els from IWooslev &: Weaver! (|1995l ) for the first second 
post bounce, bot h using the pure hadronic EoS from 
iShen et all (I1998D . Temperatures and baryon densities 
reach from below 0.1 MeV and 10 g/cm^ (6x 10~^^ fm~^) 
up to 100 MeV and 10^^ g/cm^ (0.6 fm'^). The electron 
fraction is represented on a color-scale, where low values 
of Ye — 0.1 are obtained close to the neutrinospheres 
at sub-saturation densities. For the low density do- 
main, which is entirely given by the progenitor, matter is 
isospin symmetric with an electron fraction of Yg ~ 0.5. 
Furthermore, the figures show that only a narrow mid- 
dle band of the (T, p)-diagram is covered in core-collapse 
supernova simulations of massive stars, e.g. low temper- 
atures and high densities (which corresponds to neutron 
star matter) as well as vice versa are not obtained. This 
is because the entropy per baryon does not change over 
many orders of magnitude during the early (up to 1 sec- 
ond) post bounce evolution which is illustrated in Fig. [TJ 

2.1.1. Hadronic matter 

For temperatures roughly below 0.5 MeV, the rates for 
nuclear reactions are important and the assumption of 
nuclear statistical equilibrium (NSE) cannot be applied. 
The baryon contributions to the EoS are given by nuclear 
abundances which must be evolved in time. The dynami- 
cal evolution of these abundances are c alculated using the 
nuclea r reaction network developed bv lThielemann et al.l 
(pool (and refer ences therein) which employs tabulated 
reaction rates. Due to current computational limitations, 
we employ the free nucleons, and in addition "^He and ^He 
as well as the 13 symmetric nuclei starting from ^^C to 
^^Ni. Furthermore, we use the following 3 asymmetric 
nuclei ^'^Fe, ^^Fe and ^^Fe. The details of the implemen- 
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tation of the network into AGILE-BOLTZTRAN as weh 
as the motivation for our choice of nuclei and the use in 
core-co llapse supernova models, is given in lFischer et alj 
(|2010b( ) . It is of special importance to mimic the energy 
generation from nuclear burning processes in explosion 
models, where the explosion shock travels through the 
low temperature and density envelop of the progenitor 
and the simplification of the ideal gas of Si-nuclei cannot 
be used. 

For higher temperatures, matter is assumed to be in 
NSE and the production and destruction of nuclei are in 
thermal and chemical equilibrium with respect to strong 
interactions. The abundances depend only on the ther- 
modynamic state given by the temperature T, the baryon 
density p (alternatively Hb) and the electron fraction 
Yp. For t h is reg ime, we use the hadronic EoS from 
iShen et "aTl ()1998l ) . It is based on the rclativistic mean 
field (RMF) approach for homogeneous matter and the 
Thomas-Fermi-approximation for a single average heavy 
nucleus. It has an incompressibility modulus of 281 MeV 
and an asymmetry energy of 36.9 MeV. For this stiff EoS, 
cold neutron stars with maximum mass of 2.2 Mq are ob- 



/ii, and mass rrii as follows 



tained (see the black solid line in Fig. 7(a) ). 

On top of the baryons, contributions from elec- 
trons and positrons as well as photons and coulomb 
corrections (only for no n-NSE) are added, based on 
iTimmes fc ArnettI (|1999l ). 



2.1.2. Quark matter in the hag model 

Bag models are phenomenological models which were 
originally introduced to describe quark confinement. Due 
to their simple handling an d ability to re produce hadron 
properties (see e.g. ,,Cho do s et al.iri974bl : Dcgran d et al.l 
[l97a iDetar fc Dohoghuelll983D . thev are also often ap- 
plied for bulk quark matter and phase transitions in com- 
pact star int eriors. The f i rst de scription of a bag model 
goes back to lBogolvubovl (I1968D an d was improved a few 
years later bv lChodos et al.l (|1974bl ) , known today as the 
MIT bag model. The main idea of quark bag models is 
that the true vacuum of QCD is a medium which refuses 
the penetration of quarks and confines them within a 
sphere (or a bag) to form a color neutral hadron. The 
QCD vacuum exerts a pressure on the hadron, repre- 
sented phenomenologically by the bag constant B, which 
is opposed by the motion of the interior quarks. Within 
the bag, quarks are assumed to move in asymptotic free- 
dom with vanishing or current masses. Their total pres- 
sure p*^, energy density e*^, entropy density s*^, and 
baryon number density can be written as follows 



(2) 
(3) 



where the sums run over all present quark flavors i. In the 
simple bag model quarks are treated as non-interacting 
fermions and their contributions in the above sums can 
be calculated by solving the corresponding Fermi inte- 
grals for given temperature T, quark chemical potential 



1 

Pi{mi,T,fii) = --^ J k^dkk 



dE,{k) 
dk 



X [/(fc,/i,) + /(fc, -/^,)] 



(4) 
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€i{m^,T,^i) = -^ I E^{k)k'^dk 

x[f{k,p,) + f{k,~^i,)] (5) 

Si{mi,T,iJ.,)^^ J k'^dk[-f{k,fi,)\nf{k,pi) 

-(l-/(fc,A^,))ln(l-/(fc,M.)) 

— f{k, — /ii)ln/(fc, ~pi) 

- (l-/(fc,-/i,))ln(l-/(fc,-M«))] (6) 



ni{mi,T, Pi) = 



k^dk 



27r2 ^0 

x[/(fc,MO-/(fc.-/^0]- (7) 

/(fc, ±/Xi) are the Fermi distribution functions with 
chemical potentials for particles {-\-pi) and antiparticles 



/(fc,±/i.) = 



1 



(8) 



where k is the momentum and Ei(k) = yj mf + fc^ is the 
quark Fermi energy. The number of degrees of freedom 
for each flavor gi consists of 2 spin states and 3 colors. 

Various approaches have been introduced to ex- 
tend and improve the simple bag model (see e.g. 
IDetar fc Donoghu^ Il983f l , including first order correc- 
tions for the stro ng coupling constant (see e.g. 
iFarhi fc Jaffelfl98^ . Since analytical expressions for the 
thermodynamic potentials at finite temperature and in- 
cluding as corrections can only be obtained for massless 
quarks, we calculate the pressure of massive quarks for 
the flavor i as follows 

Pi{m^,T, p^,as)=Pi[mi,T, p^,0) 

+ [Pi{0,T, Pi,as) -pi{0,T, p,,Q)] 



--Pt{m^,T,pi,0) 



7 4 250c^ 

1 TT 

60 2l7r 



^ ( It-,,- + ^ 



(9) 



The last two terms in Eq.Q are taken from analytical 
expressions in lFarhi fc Jaffj (|l984[ ). where pi(mi, T, pi, 0) 
can be calculated by numerically solving the Fermi inte- 
grals. The energy density e^, number density n^, and en- 
tropy density Si can be calculated in a similar way. We 
will apply this procedure only for the strange quarks, 
for which we choose a mass of nig = 100 MeV in ac- 
cordance with the range of nig ^ 70-130 MeV and the 
weighted average of lOStij MeV fromjAmsld (|2008l) . 
The up and down quarks have masses of several MeV 
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and can be treated as massless, that is rriu — md ~ 0. 
The remaining quark flavors are too heavy to appear in 
supernova and neutron star environments. 

The value of the bag constant is an active subject of 
research and expected to be i n the range of ^1/4 145- 
235 MeV, from hadron fitting (jPetar fc DonoghudHOSl . 
An upper limit is difficult to define. However, concern- 
ing compact stars a value of i?^/^ < 200 MeV allows the 
presence of a s mall pure quark matte r core in the star 
interior (see e.g. iSchertler et"all |2000[) . The lower limit 
for the bag constant depends on whether strange quark 
matter is considered as the ground state of nuclear mat- 
ter or not. T his absolute stability was introduced by 
iWittenI (|1984[ ) and is based on the idea that hadronic 
matter is a metastable state while strange quark matter 
has a lower energy per baryon than ^^Fe and is therefore 
the true ground state of nuclear matter. This so-called 
Witten hypothesis leads to the existence of strange stars, 
which are composed of absolutely stable strange quark 
matter. These can either represent all compact stars or 
be in co-exi stence with hadronic s tars, as was recently 
discussed bv lBauswein et afl (|2009[ ). 

In the simple bag model with m„ = = and 
nis = 100 MeV, the absolute stability of strange quark 
matter sets in for B^ ^^ < 161 MeV. For fin ite and dif- 
ferent quark masses, iFarhi fc JafFd (119841) have mapped 
out t he limiting values of B (see also IWeissenborn et al.l 
120 111 for different as)- In the present article, we will not 
assume that strange quark matter is absolutely stable. 
However, since our aim is to probe low critical densi- 
ties for the quark matter phase transition in the early 
post-bounce phase of supernovae, we chose the parame- 
ter sets of 5^/4 = 162 MeV (EOSl) and B^/^ = 165 MeV 
(E0S2) both with = and B^/'^ = 155 MeV with 
a, = 0.3 (E0S3). 

2.2. The hybrid EoS 

To probe the appearance of quark matter in supernova 
environments we implement the quark-ha dron phase 
transi tion into the pure hadronic EoS from iShen et al.l 
()1998| ) . The critical density for the onset of quark matter 
and the properties of the mixed phase depend strongly 
on the choice of global or local conservation laws. In this 
work we will apply the Gibbs approach as discussed in 
iGlendenningl (|1996f ) , where conservation laws are always 
fulfilled globally and the pressure in the mixed ph ase is a 
smooth function of the density (see Figs. 2(a) and 2(b) ). 

Another frequently used method to calculate phase 
transition is the Maxwell construction. Here only the 
baryon number is a global quantity and other conser- 
vation laws must be fulfilled locally in the quark and 
hadronic phases. As a consequence of this restriction, 
the onset of the mixed phase for the Maxwell condition 
is generally later (i.e. at higher density). Furthermore, 
for an isothermal phase transition the pressure is con- 
stant throughout the entire mixed phase, with respect to 
the baryon density. 

Despite their different behaviors, the Maxwell and the 
Gibbs approaches both represent first order phase tran- 
sitions. Gibbs falls in the class of non-congruent phase 
transitions whereas the restriction to only one globally 
conserved quantity in the case of a Maxwell construc- 
tion is congruent (for details on (non) congruent phase 



transitions, see e.g. Ilosilevskivl l2010f ). Similar to the 
pasta phases in nuclear matter at the liquid-gas phase 
transition, finite size effects from e.g Coulomb interac- 
tions and surface tension of quark or hadronic matter 
can lead to the formation of structures in the mixed 
phase, like spheres, rods, an d planes. It was shown 
by e.g. lYasutake etall (|2009[ ). that such a scenario can 
make the quark-hadron phase transition very similar to 
the one given by a Maxwell construction. However, in 
the present article we neglect any finite-size effects and 
Coulomb contributions and treat the phases in the ther- 
modynamic limit. 

Both the Maxwell and the Gibbs approaches are used 
in astrophysics and nuclear physics frequently. However, 
they represent only two extreme scenarios for equilib- 
rium cond itions at the phase tr ansition. It was recently 
shown by iHempel et al.l (j2009() that the actual variety 
of possible combinations of global and local conservation 
laws is much richer and can open up interesting new sce- 
narios for e.g. the ea rly evolution of protoneutron stars 
(jPagliaraet all [20091 ). 

Chemical potentials which correspond to a globally 
conserved quantity k, arc equal in the quark and hadron 
phases. Together with thermal and mechanical equilib- 
rium of quark and hadronic matter, the conditions for 
the co-existence region are: 



(10) 

(11) 
(12) 

The pressures p^ and p^ , in the hadronic and quark 
phases, respectively, are in principle sums over contribu- 
tions from all present particles, including electrons and 
neutrinos. As discussed in §2.1, the high temperatures 
and densities when quark matter sets in, correspond 
to conditions where neutrinos are completely trapped. 
Hence, neutrinos are in weak equilibrium with nuclear 
matt e r at given lepton frac tion (see e.g. lSteiner et al.l 
120001 : IPagliara et al.l [2009t) . together with the electrons 
and positrons. However, within the Gibbs approach the 
electron and neutrino chemical potentials are equal in the 
quark and hadronic phases and so are their contributions 
to the thermodynamic quantities such as the pressure, 
which consequently cancel in the equilibrium condition 
of Eg.JT^. In this c ontext, it was recently discussed by 
IHempel et al.l ()2009() that neutrino contrilDutions do not 
in fact have to be taken into account, neither for the 
Gibbs nor the Maxwell construction if the EoS is pro- 
vided for given temperature T, baryon density ub and 
proton fraction Yp. 

The chemical potentials for the up and down quarks 
can be obtained from the equilibrium conditions at de- 
confinement as follows 



2u + d, n ^ 2d + u, 



where 



Mil ~ 3^"' 



(13) 



(14) 



The dynamical timescales considered in supernova sim- 
ulations can be as small as milliseconds. On the 
other hand, weak reactions producing strangeness oper- 
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Fig. 2. — Hybrid EoS with a global proton fraction of Vp = 0.3. (a) shows the pressure versus density for fixed entropy per baryon of 
s = 2 kfl, for different bag constants B and ols- (b) shows pressure verses density for different entropies per baryon at fixed bag constant 
of Sl/"* = 165 MeV and ag = 0. Due to the smaller critical density, the EoS in the co-existence region is softer for higher values of s. 

Furthermore, within the Gibbs approach, the proton 
fraction Yp is a global quantity and composed of the 
charge fractions in the hadronic and the quark phases 
as follows 



ate much faster on the order of 10~^ s (from kaon de- 
cays; lAmsleii I2008D . The latter should be equilibrated 
in nuclear matter and hence be present either already 
in the hadronic ph ase in form of hyperons (see e.g. 
Ilshizuka et alj [20081 ) or be produced by a series of weak 
reactions in the quark phase, that can be schematically 
expressed by the following effective reaction 



Y^UB = (1 - x^y"^ 



H 



X 



(20) 



u + s , 



(15) 



which leads to /i^ = ^.d for the strange quark chemical 
potential. Consequently, we consider quark matter com- 
posed of the three flavors up, down, and strange. We 
did not implement hyperons or kaons in the hadronic 
phase, because our main goal is to study the effects 
from quark matter in core-collapse supcrnovae. However, 
such inclusion of strange hadrons, also in connection with 
strangeness conserving quark matter phase transitions, is 
of high interest and should be addressed in future studies. 

The mixed phase is characterized by the fraction of 
quark matter x, as follows 



where ~ np/n% (hadronic matter) and Y^ = 
(quark matter) are calculated from the positively charged 
up and negatively charged down and strange quarks via 
= (2n„ - lid - n,)/3. 

2.3. Discussion of the Hybrid EoS 



(16) 



where V'^ and are the volume fractions of the quark 
and hadronic phases, respectively. Consequently, 



<x< 1 



(hadronic phase), 
(mixed phase), 
(quark phase). 



For < X < 1 1 the baryon number density n b , the energy 
density e, and the entropy density s are calculated as a 
sum of the quark and the hadron contributions via the 
following expressions 



(17) 
(18) 
(19) 



ub = 


(1 


- X)nf - 


f Xn% 


e = 


(1 


- xy" ^ 




s — 


(1 


- x)s" ^ 


VXS^ 



Fig. 2(a) shows the pressure with respect to the baryon 
density tor different values of B and as as well as for dif- 
ferent entropies per baryon. The EoSs are calculated for 
a global proton fraction of = 0.3 at fixed entropy s = 2 
ks /baryon. The onset of the mixed phase happens for 
all parameter sets around nuclear saturation density. It 
can be identified by a softening of the EoS Q, while the 
appearance of the pure quark matter phase is accompa- 
nied by a stiffening. In the following, we will discuss the 
origin of these changes in the pressure slopes while their 
effects on the supernova dynamics will be addressed in 

§3- . 

For the simple bag model, without as corrections, 
the smaller bag constant B^/^ = 162 MeV leads to 
an earlier onset of the mixed phase and therefore also 
to a lower pressure in the co-existence region than for 
B^/^ = 165 MeV. However, as given in Eq.(I5]) and can 
be seen from Fig. [2(^ for pure quark matter the larger 
bag constant gives a lower total pressure, while the inclu- 
sion of fi rst order as correc tions stiffens the quark EoS 
(see also lAlford et al.ll2005f) . Hereby, the accompanying 
shift in the critical density to higher values can be de- 
creased again by lowering the bag constant (in our case 
B^/^ = 155 MeV) fsee lSagert eFalllMl . However, the 

^ The softness/stiffness of the EoS is measured via the adiabatic 
index 7 = 91np/91np. 
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Fig. 3. — Symmetry energy of hadronic matter from lShen et al] 
II1998I ) (solid line), 2 flavor quark matter (dashed line) and three 
flavor quark matter (dash-dotted line) v^ithin the quark bag model 
with Bi/"* = 165 MeV. 

higher pressures in the mixed and the quark phases due 
to as corrections remain. With increasing entropy, the 
critical density shifts to lower values, and furth ermor e 
the pressure in the mixed phase is reduced (see Fig. 2(b) ) 



The global proton fraction Yp also impacts the critical 
density due to the symmetry energy Esym of hadronic 
and quark matter. For nucleons, Egym is composed of 
a potential part which is caused by isospin dependent 
nucleon-nucleon interactions and a kinetic part stem- 
ming from the proton and neutron Fermi energies. In 
contrast, quark matter in the bag model is sensitive 
to its isospin state only by the Fermi energies of the 
quarks and has therefore a smaller symmetry energy 
than hadronic matter. This can be seen in Fig. [3l which 
shows Esym for hadronic matter, for two-flavor and for 
thre e-flavor quark matter with i ? ^/^ = 165 MeV (see 
also iPagliara fc Schaffner-BielichI (|2010f )). The symme- 
try energy can be defined by the difference of the energy 
per baryon of pure neutron matter, that is Yp = 0, and 
isosp in symmetric matter with Yp = 0.5 (see iLi et alj 
[20081) 



E. 



sym[nB) OC — 



E 

A 



(21) 



Tii5,yp=o.5 



Because the down and strange quarks have negative 
charge, they can replace the electrons as agents for charge 
neutrality. In Fig. [31 /3-equilibriurr|3 was assumed, which 
leads to compositions with small values of Ye or Y^ , re- 
spectively. As a consequence, we see in Fig. [3| that Esym 
is even negative for strange quark matter in contrast to 
two flav or qua rk matter and hadronic matter. 
In Fig 4(a) the charge fractions Y^ and Y^ are plot- 



ted as functions of the baryon density for the mixed phase 
as well as the regions with pure hadronic and quark mat- 
ter. The quark EoS is calculated using 

5^/4 = 165 MeV 

in /^-e quilibrium at T = 0. As described bv iGlendennind 
()1996[) . the large asymmetry energy in hadronic matter 
can be relieved by arranging the up, down, and strange 
quark fractions in the mixed phase to negative values 
of Y^ . This is also mirrored in the behavior of the 
EoS which experiences a prompt softenin g at the be- 



ginning of the mixed phase, seen in Fig. 4(a) With 



growing amount of quarks and decreasing contribution 

^ /3-equilibrium is given by the followi ng relation b etween the 
chemical potentials /le -\- /ip = fin (see e.g. IWebeiiri999l ) 
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Fig. 4. — Pressure (black solid lines) and charge fractions in the 
hadronic and quark phase, (red dash-dotted lines) and (red 
dashed lines), respectively, for matter in /3-equilibrium in graph 

(a) and at a fixed global proton fraction of 1^ = 0.3 in graph 

(b) . Both calculations are done for zero temperature and show 
the corresponding pressure-density curves, together with different 
quark volume fractions x in the mixed phase. 



from hadronic matter, the latter can approach its ener- 
getically more favorable state with Y^ ~ 0.5, while the 
charge fraction in the quark phase reaches the discussed 
small value of Y'^ ^ 0. 



As shown in Fig. 4(b) for Yp = 0.3 the EoS does not ex- 
perience a significant softening at the onset of the phase 
transition, but stays similar to the one of pure hadronic 
matter. The reason is that at 1^ = 0.3 nucleonic mat- 
ter is already close to its energetically favorable isospin 
symmetric state. Only for larger x tl^c lower pressure 
of quark matter starts to dominate and the EoS in the 
mixed phase softens. 

Due to the dependence of the hadronic EoS on Yp, the 
critical density for the phase transition is also sensitive 
to the proton fraction. The earlier onset of the mixed 
phase for smaller v alues of Yp was already present ed by 
iDrago fc Tambinil (|1999[ ) and is seen in the Figs. 5(a) 
5(b)] and [5(c)] for B^/^ = 165 MeV, B^/^ = 162 MeV and 
"B^ — 155 MeV with as = 0.3, respectively (see also 
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Fig. 5. — Phase diagrams for transitions from hadronic matter based on lShen et al.l II1998I ) (RMF, TMl) to strange quark matter with 
fixed local proton fractions of Yp = 0.1, 0.3, and 0.5 and different quark EoS parameter choices in the graphs (a)-{c). In addition, Fig. (d) 
compares a phase diagram for the conditions in supernova environments where Yp = 0.3 (strange quark matter) with heavy-ion collisions 
where the proton fraction is close to isospin symmetry and the phase transition takes place from hadronic to up and down quark matter. 
The thin lines show the onset of quark matter and the thick lines the pure quark phase. 

iFischer et aLll2010al) . 

Another important aspect for the phase transition is 
the different number of degrees of freedom of hadronic 
and quark matter. As discussed above, it is appropriate 
to assume the quark-hadron phase transition to strange 
quark matter in supernovae. The strange quarks are pro- 
duced and equihbrated by weak interactions in super- 
nova environments due to the timescalcs which are on 
the order of miUiscconds. In heavy-ion colhsions, the hot 
firebaU of nuclear matter expands within 10"^'^ seconds 
(jKolb et alJl2000f ) and thereby prevents the appearance 
of net strangeness via weak reactions. Hence, the quark- 
hadron phase transition occurs only to up- and down- 
quark matter in heavy-ion collision experiments. Fur- 
thermore, for supernovae and neutron stars we can re- 
gard phase transitions from hadronic matter to strange 
quark matter at low global proton fractions Yp <C 0.5. 
In heavy-ion collisions, the colliding nuclei are close to 



isospin symmetry with a proton fraction of ~ 0.5. In 
the case of only two flavor quark matter, the stiffening of 
the quark EoS due to the fewer degrees of freedom and 
the softness of hadronic matter makes a phase transition 
to quark matter less favorable and shi fts the critical den- 
sity to higher values. As shown in Fig. 5(d)[ a low critical 
density close to normal nuclear matter density in super- 
nova and neutron star environments is compatible with 
a high critical density (five times nuclear matter density) 
in heavy ion experiments, based on the simple bag model 
for the desc riptio n of qu ark matter. 



The Figs. 6(a) - 6(d) show the evolution of the tem- 
perature in the pure and the mixed phases along dif- 
ferent iscntropes (i.e. curves of constant entropy per 
baryon) for fixed proton f ractio n and three-flavor quark 



matter. Hereby, the Figs. |6(a) 

agrams for B^^^ 



6(c) give the phase di- 



165 MeV for different values of Y„ 
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Fig. 6. — Temperature evolution in the mixed phase for the different entropies per baryon s = 1,2, 3,4 kg. The graphs (a) - (c) show 
calculations for global proton fractions of Yp = 0.1, Yp = 0.3 and Yp = 0.5 with B^/* = 165 MeV. The red so lid lin es show the onset of 
the mixed phase (thin lines) and the beginning of the pure quark phase (thick lines) and are the same as in Fig. |5(a)| For comparison, the 
graph (d) shows the temperature evolution in the mixed phase for Yp = 0.3 and B^/'* = 200 MeV. 

Dynamical simulations (e.g. in astrophysics and heavy- 
ion collisions) that use microscopic EoSs raise the ques- 
tion about convexity, in particular if phase transitions 
are included. In the context of phase transitions, one 
very often discusses the convexity of the thermodynamic 
pote ntial (in our ca,se the Helmholtz free-energy, see 
e.g., iMiiller &: Serot I [19951 ): If the thermodynamic po- 
tential is a convex function of its arguments, the sys- 
tem is thermodynamically stable. Due to the construc- 
tion of the quark-hadron mixed phase and the stabil- 
ity of each phase separately, our hybrid EoS is stable 
and thus convex. However, another definition of con- 
vexity can be found in the literature in the context of 
the hydrodynamic Riemann problem, which refers to the 
curvature of the pressure as a function of density along 
isentropes. This property of the EoS is particularly rel- 
evant for shock waves traversing the coex istence region 
of a f irst order phase tra nsition (see e.g. iRischke et alJ 
119901 : iBugaev et aI1ll989| ) as unusual wave patterns and 
effects like wave s plitting can o ccur. Within this defini- 
tion of convexity. iHeuze et al.l (j2008i ) explored whether 
finite volume methods give a reasonable description of 



while the Fig. |6(d)] shows the temperature evolution for 
B^/'^ = 200Mey and = 0.3. ISteiner et al.1 (l2000l) and 
INakazato et al.l (j2010bD discuss such a temperature evo- 
lution along isentropes for B^^^ ~ 200 MeV at a fixed 
lepton fraction Y^. The authors find a temperature de- 
crease in the mixed phase due to the larger number of 
degrees of freedom and the different heat capacity for 

relativistic quarks. 

However . iDrago fc Tambinl (|1999f ) showed that the be- 
havior of T in the mixed phase depends on the critical 
density and the temperature at the onset and at the end 
point of the mixed phase, which is set by Yp. The tem- 
perature behavior base d on the quar k bag model can 



be s een from the Figs. 6(a) and 6(b) in comparison to 



6(c) and |6(d)| For early phase transitions, the temper 



ature curves rise in the beginning of the mixed phase, 
while for high critical densities the temperature in the 
hadronic phase has already reached large values and the 
onset of quark matter causes a prompt temperature de- 
crease along the isentropes. These results are in qualita- 
tive ag reer nent with the pre vious studies of lSteiner et al.l 
(|2000l) and lNakazato et al\ (|2010b[ ). 
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IShen et al.l ifTOQSl 'l (TM1,RMF) and hybrid stars with different bag 
constants using the Gibbs construction for the phase transition. The 
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(b) Neutron star density profiles (thick blue linos) including the 
quark volume fractions (thin red lines). 



Fig. 7. — Different mass-radius relations, comparing different bag constants, in graph (a) and neutron star density profiles, for two selected 
bag constants, in graph (b). The stiffening of the hylarid EoS when corrections from the strong coupling constant as arc included, leads to 
higher maximum masses of the hybrid star branches. Small values of B^/^ < 170 MeV as well as large bag constants of B^/* > 190 MeV can 
reach the mass limit of 1.44 Mq (horizontal solid line in graph (a)). Additionally, mass limi ts of 1.65 IsAq (PSR J1903+0327. iFVeire &: Wexl 
I20ia horizontal dashed line) and 1.97 ± 0.04 Mq (PSR J1614-2230, IDemorest et al.ll20TOl . horizontal dash-dotted line) are also shown in 
graph (a). The pure quark phase in the maximum mass configuration is large for small bag constants, while for B^/* = 200 MeV quark 
matter is present only in form of a mixed phase. 

a Rieniann problem with a non-convex EoS. After per- 
forming sueh a test we can conclude that AGILE belongs 
to the class of hydrodynamics codes, which show a sat- 
isfactory agreement with the analytical solutions. 



(IDemorest et alJ l201Clf ). respectively. As shown in 



2.4. Masses and Radii of Hybrid Stars 



Fig. 7(a) shows hybr id star mas s -radiu s relations for 
the hadron EoS from iShen et all (|1998[ ) and the bag 
model for different values of B and as- For compari- 
son, we indicate the mass of the Hulse-Taylor pulsar of 
AI ~ 1.44 M© 0. For simple quark bag models where 
the phase transition is constructed applying the Gibbs 
conditions, only the small values of B^^^ < 170 MeV as 
well as large bag constants of B^^'^ > 190 M eV can reach 
this limit. However, as shown in Fig. |7(b)[ hybrid stars 
constructed using B^/^ = 200 MeV, do not contain pure 
quark matter but only the mixed phase with quark mat- 
ter fractions of x ^ 0.5. In contrast, for the small bag 
constant i?^/^ = 165 MeV the hybrid star at maximum 
mass is almost completely composed of pure quark mat- 
ter, with only a thin hadronic crust with a thickness of up 
to one kilometer. As will be discussed in the following, 
the small value of the critical density and the resulting 
early phase transition to quark matter, together with the 
stiffening of the EoS in pure quark matter, lead to very 
specific dynamical consequences for the core-collapse su- 
pernova which may not be reached for large values of B 
due to the high critical density and absence of the stiff- 
ening of the hybrid EoS. 

Recent mass limits for physical EoSs come from the 
millisecond pulsars J1903-K0327 and J1614-2230 which 
were ann ounced to have m asses of M = 1.67 ± 
0.01 Mq (|Freire fc WexllMol) and M = 1.97 ± 0.04 Mq 

^ The Hulse-Taylor pulsar is the highest precisely known mass of 
compact stars, at present. It represents therefore a minimal limit 
for masses which should be reachable with realistic nuclear EoSs. 



Fig. 7(a) within the simple bag model such masses can 
only be reached for bag constants of B^^^ > 200 MeV. 
Further corrections from the strong interaction coupling 
constant can stiffen the normal bag EoS and therefore 
lea d to higher maximum masses of c ompact stars (see 
e.g.lSchertler et al.ll2000l:IAlford et al.li2005;,Sa gert et alJ 
120101 : iWeissenborn et al. I I2011D . It can be seen from 
Fig. 7(a) that the chosen parameter set B^/^ = 155 MeV 
and as = 0.3 gives a maximum mass of hybrid stars of 
-^j^max ~ 1-67 Mq, which w ould be co mpatible with the 
observation of J1903-H0327 (iFreire fc Wex.,201(j ) . The re- 



cent mass measurement bv lDemorest et al.l (|2010[ ) points 
to an even higher neutron star m ass of 1.97 ± . 04 M© . 
However, as discussed i n e.g. iKurkela et al.l (|2010( i: 
lAlford et all (f20?)5l [200l); lOzel et al.l (|20100 . such high 
values exclude neither the presence of quark matter in 
compact star interiors nor low critical densities. The 
stiffening effect of strong interactions on the quark EoS 
can be modeled in our approach by e.g. according choices 
of larger values of 0^5 and will be tested in the future. 

2.5. Neutrino interaction treatment in quark matter 

Based on the parameters chosen for the bag model as 
discussed in the previous subsections, quark matter oc- 
curs at nuclear matter densities. At such high densities, 
the neutrino mean free paths are extremely small and 
neutrinos are trapped. Hence, we always assume chemi- 
cal, weak and thermal equilibrium when determining the 
neutrino chemical potentials in quark matter. 

3. SIMULATION RESULTS 

In this section, we report on results from core-collapse 
supernova simulations in spherical symmetry where we 
apply the quark-hadron hybrid EoSs introduced in §2.3. 
The quark-hadron phase transition was found to oc- 
cur during the early post bounce phase within the first 
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Fig. 8. — Evolution of selected properties of the innermost mass trajectories in core-collapse supernova simulations of massive progenitor 
stars at two examples. The red dashed lines in the density graphs illustrate nuclear saturation density and the solid blue lines show the 
central data. 



500 ms for all cases considered in this investigation. This 
is the central property of the quark matter description 
applied here due to the selective choice of parameters, 
which lead to the onset of the quark-hadron phase tran- 
sition close to nuclear saturation density. Consequences 
with respect to the dynamical evolution will be further 
discussed below. The simulations are launched from low 
and interme diate mass Fe-core pro genitors of 10.8, 13 and 
15 Mq from lWooslev et all (|2002[ ). First results from the 
10.8 and the 15 s i mulat ions have already been dis - 
cussed in lSagert et all (|2009( ) and lFischer et alJ ()2010af ). 
Here, we extend their analysis and present a deeper inves- 
tigation with respect to the appearance of quark matter 
in PNS interiors and discuss consequences with respect 
to the neutrino signal emitted. 

3.1. The early post bounce phase 

The early post bounce phase of core-collapse super- 
novae of massive Fe-core progenitors is determined by an 
extended mass accretion period, during which the cen- 
tral density and temperature increase above several times 
nuclear density and above several tens of McV. This be- 
havior is shown in Fig. [8] for the example of the low- 
mass 10.8 M0 and the high-mass 40 M© prog enitor mod- 
els, bo th using the (pure hadronic) EoS from iShen et all 
()1998() . Illustrated is the time evolution of selected mass 
trajectories, which belong to the PNS interiors, during 
the Fe-core collapse (see the radial evolution), bounce 
and the first 400 ms post bounce. Both models reach sev- 
eral times nuclear saturation density (red dashed lines in 
the density graphs) at bounce (i.e. at t = 0). They pro- 
ceed in a similar fashion during the post bounce compres- 
sion after the expanding bounce shock has come to halt. 
The central densities increase slowly above nuclear satu- 
ration density and reach 3.36 x 10^ g/cw? (0.201 fm~^) 
for the 10.8 M© progenitor model and 3.99 x 10"'^"' g/cm^ 
(0.238 fm""^) for the 40 M© progenitor model within the 
first 500 ms after bounce. Furthermore, the central elec- 



tron fraction reaches low values, between Ye — 0.28 (at 
the center) and Kg — 0.1 near the neutrinospheres at 
intermediate densities around 10^^ — 10"'^'^ g/cvo?. This 
highly deleptonized region originates from the launch of 
the deleptonization burst at about 5 ms post bounce 
(dominated by Ve)- Differences between these two pro- 
genitor models occur in the temperatures obtained dur- 
ing the first 400 ms post bounce evolution. The 10.8 M© 
progenitor model in Fig. |8] (a) reaches temperatures of 
T = 15 - 35 MeV and the 40 M© model in Fig. [5] (b) 
between 20 — 60 McV. This effect is due to the more 
compact PNS obtained at bounce for the 40 M© progen- 
itor model and the larger mass accretion rate at the PNS 
surface. 

3.2. The quark-hadron phase transition 

We apply the hybrid EoS introduced in § 2.3 to core- 
collapse supernova simulations of massive Fe-core pro- 
genitor stars in spherical symmetry, where no explosions 
are obtained using standard hadronic EoSs from e.g. 
iLattimer k, Swestvl ()1991| ) and lShen et~all (|1998( ). In the 
following paragraphs, we illustrate the dynamical evolu- 
tion obtained during the quark-hadron phase transition. 
We select as standard reference model, the 10.8 M© pro- 
genitor model applying E0S2. 

For this choice of parameters, quark matter appears 
at densities close to nuclear saturation density, where 
the critical density depends strongly on the temperature 
and the electron fraction obtained (see the discussion in 
§ 2.3). The conditions permit quark matter already at 
core bounce, as illustrated in Fig. [9] (left panel). The 
corresponding hydrodynamic properties obtained are il- 
lustrated in Fig. 1101 Quark matter appears at the cen- 
ter (< 0.3 M©), at a mass fraction of about x=0. 01-0.1 
(see the radial bounce profiles in Fig. [S] (left panel) and 
Fig. [TU] for baryon density, electron fraction and temper- 
ature). The figures also show that down-quarks are the 
most abundant species, followed by strange-quarks and 
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Fig. 9. — Evolution of the mass fractions of nuclcons and the different quark flavors at bounce (left panel) and at the two selected post 
bounce times 100 ms (middle panel) and 400 ms (right panel) for the 10.8 Mq reference model. 
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Fig. 10. — Radial profiles of selected hydrodynamic quantities 
during the appearance of quark matter for the 10.8 Mq reference 
model at bounce (solid lines) and at the two selected post bounce 
times 100 ms (dashed lines) and 400 ms (dash-dotted lines). The 
horizontal dotted line in the entropy graph illustrates the critical 
adiabatic index of 7 = 4/3. 

up-quarks, which can be understood as a consequence of 
charge neutrahty and the low electron fraction at the cen- 
ter. This hierarchy remains during the later post bounce 
evolution. The small central fraction of quark matter re- 
duces during the early post bounce expansion (see Fig. |9] 
(middle panel) and the radial profiles at bounce and at 
100 ms post bounce in Fig. llOp. Only after the dynamic 
bounce shock has stalled and the subsequent evolution 
is determined by mass accretion onto the PNS and heat- 
ing behind the standing accretion shock (SAS), does the 
central density increase and the central electron fraction 
decrease. This in turn causes the quark matter volume 



fraction to rise (sec Fig. [5] (right panel) and Fig. [TU] at 
400 ms post bounce) . The timcscalc for the quark matter 
fraction increase is given by the central density and tem- 
perature increase as well as by the electron fraction de- 
crease. Hence it is determined by the mass accretion rate, 
which in turn depends on the progenitor model and the 
hadronic EoS (before the EoS is dominated by the quark 
contributions). For the massive Fe-core pro genitors un- 
der in vestigation using the stiff EoS from IShen et all 
(|1998[) . the timescale is on the order of 100 ms, where 
a central density of 5.2 x 10^^ g/cm'^ and a temperature 
of 15 MeV are obtained (see Fig. (TU]). During the post 
bounce evolution, the adiabatic index 7 decreases slowly 
below the critical value of 4/3 on the same timescale, 
illustrated in Fig. [TOl Furthermore, the appearance of 
quark matter at the PNS interior proceeds adiabatically 
as shown via the entropy per baryon profiles in Fig. [TUl 
The PNS configuration is gravitationally stable as long 
as the maximum mass (given by the hybrid EoS) is not 
reached. 

A quark matter volume fraction between and 1 trans- 
lates to matter in the quark-hadron mixed phase, for 
which the adiabatic index is reduced (see §2.3). This be- 
havior is illustrated in Fig.[TT](c) at selected post bounce 
times, after the PNS exceeds its maximum stable mass 
and starts to collapse. The last stable configuration (it 
is not illustrated for simplicity) corresponds to the post 
bounce time of 428.5357 ms. The consequent soften- 
ing of the EoS for matter in the mixed phase leads to 
a dramatic dynamical evolution. The reduced adiabatic 
index causes the mixed phase of the PNS to contract, 
as shown via the radial velocity profiles in Fig. [TT] (a) 
during the collapse at three selected post bounce times. 
During the contraction, the central density and temper- 
ature increase (see Fig. [TT] (b) and (e)). Furthermore, 
weak-equilibrium is established at a lower value of the 
electron fraction due to the changed thermodynamic con- 
ditions (see Fig. [TT] (d)). The contraction accelerates 
and proceeds into an adiabatic collapse, during which 
the central entropy per baryon stays constant as shown 
in Fig. [TT] (f). The central part of the PNS coUapscs 
sub-sonically and the outer part collapses supcr-sonically 
(similar to the Fc-core collapse). The timescale reduces 
from 100 ms to milliseconds. The compression results in 
significantly higher densities and temperatures, which in 
turn favors quark matter over hadronic matter. In this 
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Fig. 11. — PNS collapse of the 10.8 Mq progenitor model due to 
the presence of quark matter, at three different post bounce times 
(solid lines: 429.2 ms, dashed lines: 429.6 ms, dash-dotted lines: 
429.7 ms) after the PNS has become gravitationally unstable. The 
same configuration as Fig. [9] but in addition graphs (g) and (h) 
show the mass fractions of nucleons (thick red lines) and quarks 
(thin blue lines) and the mean energies for (thin lines) and for 
(/x/T)-(anti)neutrinos (thick lines), respectively. 



sense, the quark-hadron phase transition which started 
slowly on timescale of 100 ms proceeds rapidly during the 
PNS collapse until the pure quark matter phase at the 
PNS center grows and a quark core of about 1.3 Mq 
is formed (see the composition in Fig. [TT] (g)). The 
quark fractions of up-, down- and strange-quarks rise 
equally whereas in the pure quark-matter phase down- 
quarks and strange-quarks are favored over up-quarks 
(see Fig. [S]). This is due to the fact that the chemical 
potentials are equal for s-quarks and d-quarks (see § 2.3) 
and matter is neutron-rich. The different symmetry en- 
ergy for quark matter, together with the higher densities 
and temperatures obtained during the PNS collapse, re- 
sult in a different weak-equilibrium where Y^. is generally 
lower. The lower Ye obtained (see Fig. [TT] (d)) between 
0.8 — 1.3 Mq reduces the pressure of the degenerate elec- 
tron gas and hence softens the EoS additionally, which 
in turn supports the collapse. In addition, the mean neu- 



trino energies are also shifted to higher values, where the 
electron flavor neutrinos are most sensitive to density 
variations. As shown in Fig. [TT] (h), the central mean 
energies of the electron neutrinos increase from about 
250 MeV before the PNS collapse to about 300 MeV. 
The largest increase of the mean electron neutrino en- 
ergy was found for the infalling material where density 
and temperature increase most, from 10"'^'* g/cm'^ to 10^^ 
g/cm^ and from 15 MeV to 45 MeV (see Fig[n] (b) and 
(e)). There, the mean neutrino energy increases from 
about 75 MeV to about 200 MeV for i/e and to about 
150 MeV for the (f^/^, i/^/^). A less pronounced increase 
of the mean neutrino energies can be found for the (^/r)- 
(anti)neutrinos, which rise from about 50 MeV to about 
60 MeV at the center. The mean free paths of all neu- 
trinos are extremely small, on the order of few 10 — 1000 
cm. Hence the neutrinos are highly trapped and cannot 
escape on the timescales on the order of milliseconds up 
to hundreds of milliseconds that are found during the 
post bounce PNS evolution. 

As discussed in §2.3, the adiabatic index rises again 
in the pure quark phase. This behavior is illustrated in 
Fig. [TT] (c) during the PNS collapse. The obtained quark 
core and the consequent stiffening of the EoS for mat- 
ter in the pure quark phase at the PNS interior, causes 
the PNS collapse to halt. A strong hydrodynamic shock 
forms, which can be identified in the radial velocity pro- 
file at t = 429.6 ms post bounce in Fig. [TT] (a) and the 
entropy increase in Fig. [TT] (f ). The shock wave forms at 
about 1.25 Mq. The system does not overshoot its hy- 
drostatic equilibrium configuration, because quark stars 
(or hybrid stars with extended strange quark cores) are 
energetically self bound objects, the total internal energy 
becomes larger than the gravitational binding energy, at 
the densities co nsidered here (for a detailed discussion, 
see e.g. § 18 in IWeberi[T999l) . Hence, the scenario can- 
not be considered a second bounce. Furthermore, when 
the shock reaches the sonic point it still remains a pure 
accretion front with no matter outfiow, as shown in the 



Figs. [TT] and |l2(a) 

The shock forms due to the stiffening of the EoS in the 
pure quark phase and due to the supersonically infalling 
outer part, which does not know about the halted col- 
lapse at the center. Information about the central stiff- 
ening cannot propagate outward across the sonic point. 
Initially, the second shock wave appears as a standing 
accretion front with no matter outflow (see Fig. [TT] (a)). 
The shock evolution is given by the balance of ram- 
pressure ahead of the shock from the infalling nucleons 
and the thermal pressure from the dissociated quarks 
behind the shock. The shock position determines the 
phase boundary between the quark-hadron mixed and 
the pure hadronic phases, due to the large density and 
temperature jumps at the shock front. The dissociated 
nuclear matter, which accumulates onto the PNS sur- 
face, has been shock heated at the SAS previously. As 
this dissociated material crosses through the second ac- 
cretion shock, it is converted from hadronic matter to 
quark matter. For all our models, the second accretion 
shock was found to propagate outward towa rds the PNS 
surface (see the radial velocity profiles in Fig. 12(a) (a) at 
the example of the 10.8 Mq progenitor reference model). 
Thereby, the outwards propagating shock wave remains a 
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(a) The initial explosion phase. 




Fig. 12. — Radial profiles of selected hydrodynamic variables during the explosion of the 10.8 M0 reference model due to the quark 
hadron phase transition, at selected post bounce times. 
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pure accretion front with no matter outflow, as shown in 
Fig. [TT](a) betwee n 429.6 ms and 429.7 ms post bounce 
and in Fig. Il2(a)| at 429.68 ms and 429.83 ms. 

The density and temperature jumps at the shock in- 
crease, as long as it remains a pure accretion front. Mass 
from the outer gravitationally unstable part of the PNS 
continues to fall onto the shock front, increasing the infall 
velocities at the sh ock (see the velocities in the middle 
panel of Fig. 12(a) I. The material accumulates onto the 
central quark core. The increasing thermal pressure be- 
hind the standing shock front drives the shock slowly out- 
wards, supported by the increasing heating right behind 
the standing accretion shock illustrated in Fig. 13(a~ 



In the presence of the largely enhanced i?e-luminosity, 
Pe-absorption at protons dominates over i/e-absorption 
at neutrons, even though the conditions are generally 
neutron-rich. The sharp rise of the heating rates just 
behind the shock corresponds to the increasing Yg illus- 
trated in the Fig. 12(a) Behind the heating region, small 
(in comparison to t he tiea ting at the shock) cooling rates 
dominate (see Fig. 13(a)). 

As the density of tlie infalling material decreases it 
becomes too low for hadronic matter to be converted 
into quark matter. At the phase boundary between the 
quark-hadron mixed phase and the pure hadronic phase, 
the density decreases from about 10^^ g/cm^ to about 
10^'' g/c\r? (see the lower panel of Fig. 12(a)). Eventu- 
ally, the shock wave detaches from the surface of the PNS 
and accelerates. Positive matter velocities are obtained, 
where where the accretion front turns into a dynamic 
shock wave. Because the density at the PNS surface de- 
creases also at larger radii, from 10^^ g/cvp? at the phase 
boundary between the mixed and the hadronic phases to 
10^° g/cm^ at the low den sity en velope, the shock contin- 
ues to accelerate (see Fig. 12(b) |. Velocities on the order 
of 1 — 1.5 X 10'^ km/s are obtained during the shock pas- 
sage across the low density envelope. The scenario is 
again different from the early shock propagation of the 
bounce shock, where the shock suffers immediately after 
its formation from the dissociation of infalling heavy nu- 
clei which causes an energy deficit of about 8 MeV per 
baryon. Here, the infalling matter is already dissociated 
and composed of only free nucleons and light nuclei. In 
this sense, the shock does not loose energy during the ini- 
tial propagation. Even more, during the initial shock ex- 
pansion neutrino hea ting st ill deposits energy behi nd the 



shock wave (s ee Fig. 13(a) and compare wit h Fig. 12(a) 
as well as Fig. 13(b) and compare with Fig. 12(b) ) 



The expanding dynamic shock wave finally merges 



with the SAS from the Fe-core bounce (see Fig. 12(b) 
at t = 430.24 ms post bounce), which remained unat 
fected from the events inside the PNS at a radius of 
about 80 km. Maximum matter outflow velocities of 
1.6 X 10^ km/s are obtained. The initial shock expan- 
sion slows down when the dynamic shock wave reaches 
infalling heavy nuclei from the outer layers after merg- 
ing with the standing accretion shock from core bounce. 
The matter velocities decrease continuously during the 
later evolution to typically 4 — 6 x 10"* km/s (depend- 
ing on the progenitor model). However, for all models 
under investigation the expanding shock wave was never 
found to stall again at later times. Hence, the quark- 
hadron phase transition triggers not only the formation 



of a strong additional shock wave but also launches an 
explosion in core-collapse supernova simulations where 
otherwise no explosions could be obtained. 

3.3. The neutrino observables 

First, we describe the standard neutrino signal ob- 
tained from a failed core-collapse supernova explosion of 
10.8 M p) using the pure hadronic EoS from iShen et "aD 
(|T998h . Below, we will compare these neutrino spectra 
with the spectra obtained taking QCD degrees of free- 
dom into account. 

The evolution of the neutrino luminosities and root- 
mean-square energies are shown in Fig. [14] (a) for the 
10.8 Mq progenitor model using the standard hadronic 
description of nuclear matter, measured at a distance of 
500 km in the co- moving reference frame (for a definitio n 
of the neutrino observables, see iFischer et aLll2010bD . 
During the deleptonization burst at bounce, the elec- 
tron neutrino luminosity rises up to several 10^^ erg/s 
on a timescale of 5 — 20 ms post bounce. 9^ are pro- 
duced only after bounce due to the decreased degener- 
acy which allows, next to electrons, also for positrons 
and the corresponding charged current reactions as well 
as pair-processes. The same holds for (/x/T)-neutrinos, 
which are produced only via pair-processes after bounce. 
The early (on timescale of 100 ms post bounce) evo- 
lution of the neutrino spectra is given by the balance 
of mass accretion and diffusion at the neutrinospheres. 
After the deleptonization burst has been launched, the 
electron-flavor neutrino luminosities increase slightly un- 
til about 150 ms post bounce, due to the moderately 
large mass accretion rate. The (/x/r)-neutrino luminosi- 
ties decrease constantly with time because their spec- 
tra are given by diffusion rather than by accretion. The 
electron flavor luminosities follow the same decreasing 
behavior after the mass accretion rate at the neutri- 
nospheres drops below a certain threshold (depending 
on the progenitor model). At 400 ms post bounce, the 
electron- flavor ((/i/T)-neutrino) luminosities reach values 
of about 5 x 10^^ (3 x 10^^) erg/s. The mean energies of 
the electron-flavor neutrinos increase continuously with 
respect to time after bounce, from 10 (15) MeV after the 
deleptonization burst has been launched to 14 (17) MeV 
at 400 ms post bounce. The mean energies of the (/i, r)- 
(anti)neutrinos decrease slightly with time after bounce 
and reach about 19 MeV at 400 ms (see Fig. [T4l(a)). 

The evolution of the neutrino observables from the sim- 
ulations applying the hybrid EoS introduced in §2.3, is 
shown in Fig.[T3](b) for the 10.8 Mq progenitor reference 
model. Compared to the standard hadronic scenario, we 
find very similar spectra for post bounce times before 
the PNS collapse (see Fig. [14] (b)). Differences occur in 
slightly higher neutrino luminosities, which are due to 
the more compact PNS where quark matter is present. 

The propagation of the second shock wave across the 
neutrinospheres releases an additional burst of neutri- 
nos, where neutrinos of all flavors carry away energy on 
the order of several 10^^ erg/s on a short timescale on 
the order of a few milliseconds. The scenario is sim- 
ilar to the bounce shock propagation across the neu- 
trinospheres where the deleptonization burst is released 
which, however, appears only in Vg. In order to under- 
stand the release and the nature of the second burst as 
shown in Fig. [14] (b) at about 430 ms after bounce, we 
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Fig. 13. — Total net heating/cooling rates for Ve, (solid lines) and Ve (dashes lines) at selected post bounce times. Negative (positive) 
values correspond to a net cooling (heating). The vertical dash-dotted lines show the position of the expanding explosion shock. 



must consider what happens when the expanding sec- 
ond shock wave passes through the formerly dissociated 
and now again shock heated nucleons. Material at the 
PNS surface has an electron fraction of Ye — 0.1. This 
low stems from the del cptoniz ation burst shortly af- 
ter core bounce (see Fig. 12(a) I. The heating b y the 



sho ck lifts the electron degeneracy (compare Figs. 12(a) 
and 12(b) ). It allows for the creation of electron-positron 



pairs toUowed by positron captures on neutrons that in- 
crease the electron fraction. This can be clearly identi- 
fied via the reduced degeneracy of the infalling material 
ahead of the shock during the initial shock expansion, 
as long as the shock front remains a standing accretion 
shock (see Fig. [12] (a)). The situation changes when the 
shock accelera tes at the PNS s urface to positive veloci- 



ties (see Figs. 12(a) and 12(b)|) . The degeneracy of the 
infalling material reduces significantly at the shock front 
(see Fig. [15] (b)). On the other hand, the degeneracy of 
the matter that accumulates behind the shock increases. 
Hence, the electron fraction of the infalling material de- 
creases while the electron fraction behind the shock in- 



creases (see Fig. 12(a) I. Note that the same behavior for 



the electron degeneracy discussed here and shown in the 
Figs. ITS] (a) and (b) (left panel), holds for the charged 
chemical potential — /Xp, which is shown in the Figs. llSl 
(a) and (b) (right panel). The charged chemical poten- 
tial decreases ahead of the expanding shock front and 
increases behind it. 

The shock passage across the neutrinospheres (see the 
vertical lines in Fig. [15]) releases the second neutrino 
burst. The increasing temperature enhances the pro- 
duction of electron-positron pairs, where most of the 
positrons are captured by neutrons producing a burst 
of electron antineutrinos. Some positrons interact via 
pair process (7) from Table [TJ contributing to all neu- 
trino species. This explains why the second neutrino 
burst is dominated by while and (i-'^/n ^'^i.r) have 
similar luminosities (see Fig. [TJ] (b)). The electron 
(anti)neutrino luminosities increases from 0.2611 x 10^'^ 
(0.2708 X 10^3) gj.g/g to 0.55 x lO^^^ (^2.2 x lO^^) gj-g/g 
a short timescale on the order of milliseconds. The same 
increase in luminosity holds for the (/x/r)-neutrinos, 
which reach about 1.5 x 10^"^ erg/s. The sharp peak of the 
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(b) Simulation using the hybrid EOS2. 
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Fig. 14. — Evolution of the neutrino luminosities and mean neutr i no en ergies for the 10.8 Mq progenitor model with respect to time 
after bounce, comparing the standard hadronic EoS from lShen et al] {1993) and the hybrid EOS2. The data are measured at a distance of 
500 km in the co-moving reference frame. 



neutrino luminosities is accompanied by a sharp rise of 
the mean neutrino energies. The electron (anti)neutrino 
mean energies increase from 16.5 (18.4) MeV to 29.38 
(32.2) MeV, and the mean (/i/r )-ncutrino energies rise 
from about 18.5 MeV to 53.7 MeV on the same short 
timescale on the order of milliseconds. 

The later evolution of the PNS surface is shown in 
Fig. [16] via the radial profiles of p, Yg, velocity and T 
at selected post bounce times during the ongoing explo- 
sion phase. Neutrino cooling after the sh ock pr opaga- 



13(b) I 



leads 
ex- 



tion across the neutrinospheres (see Fig, 

to the establishment of matter infall between the 
panding explosion shock and the PNS surface already 
at about 430.24 ms after bounce in Fig. |12(b)| between 
10-20 km. As matter continues to fall onto the PNS 
surface with even supersonic velocities, an additional ac- 
cretion shock forms at the PNS surface. This additional 
accretion shock on top of the PNS surface expands due 
to neutrino heating and falls back due to cooling (sec 
ther efore t he large cooling rates between 7-10 km in 



Fig. 13(b) ) on timcscalcs on the order of 10 ms, during 



which the neutrinospheres expand and contract corre- 
spondingly. The same phenomena apply as discussed in 
the paragraph above. The shock propagation towards the 
neutrinospheres increases the temperature at the neutri- 
nospheres, where the enhanced opacities cause the prop- 
agation of the neutrinospheres to lower densities. Neu- 
trino cooling, as well as the obtained matter expansion at 



the neutrinospheres, shifts the neutrinospheres to higher 
densities. The consequently increasing and decreasing 
mass accretion rate is reflected in the electron flavor neu- 
trino luminosities and mean energies, shown in the insets 
in Fig. [14] (b) after the second neutrino burst between 
432-500 ms. 

The additional standing accretion shock, propagating 
back and forward at the PNS surface, settles down to 
a quasi-stationary state at about 550 ms post bounce. 
After that, the evolution of the neutrinospheres is deter- 
mined by mass accretion on the order of few 0.1 Mq/s 
at a timescale on the order of 100 ms up to seconds. The 
neutrino luminosities decrease to values below 10^^ erg/s 
within the first 500 ms after the onset of the explosion. 
The mean neutrino energies follow the same behavior. 
The mean electron (anti)ncutrino energies decrease from 
about 18 (20) MeV at 550 ms post bounce to 12 (18) at 
1 second post bounce (see Fig. [14] (b)). The mean (^/r) 
neutrino energies decrease only slightly from 22 MeV at 
550 ms post bounce to about 19 MeV at 1 second post 
bounce. Unaffected by the behavior of the standing ac- 
cretion shock at the PNS surface, the explosion shock 
continues to expand to larger radii with matter veloci- 
ties on the order of several 10^ km/s (see Fig. [TB]). 

4. DISCUSSION 

In this section, we will compare the results obtained us- 
ing the 10.8 Mq reference model with the remaining sim- 
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(a) The early shock propagation. 
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Fig. 15. — Radial profiles for the degeneracy parameter rj = fie/T 
and charged c hemica l po tential ^„ — for the selected post bounce 
times of Figs. |12(a)| and |12(b)| 

ulations that have been performed where different pro- 
genitor models as well as different quark EoS parameters 
are used. Further below, we will discuss possible appli- 
cations with respect to the PNS structure and the for- 
mation of magnetars as well as with respect to possible 
nucleosynthesis investigations . 



TABLE 2 

Summary of the models under investigation. 



Prog. 


EoS 






rpa 




[Mq] 




[ms] 


[W* g/cm3] 


[MeV] 




10.8 


EOSl 


240 


6.607 


13.14 


0.234 


10.8 


EOS2 


428 


6.457 


14.82 


0.237 


13 


EOSl 


235 


6.493 


13.32 


0.235 


13 


EOS2 


362 


7.228 


16.38 


0.191 


15 


EOSl 


172 


7.523 


17.15 


0.170 


15 


EOS2 


275 


7.586 


16.25 


0.187 


15 


EOS3 


308 


5.511 


17.67 


0.197 



" Selected central properties of the PNSs at the onset of collapse 



4.1. Explosions in spherical symmetry 

In addition to the low-mass reference model of 
10.8 M0, we apply the hybrid EoSs introduced in § 2.3 
to core-collapse supernova simulations of the intermedi- 
ate mass Fe-co r e prog enitor stars of 13 and 15 Mq from 



IWooslev et al.l (j2002 ). The results are summarized m 
Table [H We list the characteristic properties of the sim- 
ulations for the different progenitor models and the hy- 
brid EoSs with the two different values of the bag con- 
stant. The values tpb are the post bounce times when 
the PNSs become gravitationally unstable and start to 



collapse. The corresponding central conditions, i.e. den- 
sity pc temperature Tc and electron fraction Y^, are also 
listec0. The models using the larger bag constant, which 
corresponds to a higher critical density, reach typically 
lower quark volume fractions at equal evolutionary states 
post bounce. In order to reach similar central densities, 
the PNSs have to accrete more mass which leads to a 
longer accretion time post bounce. This effect is com- 
pensated by the lower maximum stable mass of the con- 
figurations using the larger bag constant. Hence, the 
central densities required for the PNSs to become grav- 
itationally unstable and collapse are rather similar for 
the same progenitor model. On the other hand, due to 
the longer mass accretion period for the models using 
the larger bag constant on the order of 100 — 200 ms. 
the PNSs are more compact before collapsing. In ad- 
dition, the central temperatures obtained are higher by 
about 2 — 3 MeV. This in turn favors quark matte r ove r 



hadr onic matter (see the phase diagrams in Fig. 4(a) 
5(c) in § 2.3), where weak-equilibrium is established at a 



lower value of the electron fraction (see Table [2|). 

All models listed in Table [2] evolve in a similar fashion 
during the quark-hadron phase transition. The explo- 
sions are obtained due to the formation of the strong 
second hydrodynamic shock front inside the PNSs. The 
shocks accelerate at the PNS surfaces, i.e. the shock 
breakout, which triggers the explosions where otherwise 
no explosions could have been obtained in spherical sym- 
metry. The explosion energy estimates E^xpi and ap- 
proximate neutron star masses Mns (for a definition of 
the ma ss cut and the explosion energy, see iFischer et alj 
l2010bl following a suggestion by S. Bruenn), are listed 
in Table [3l Moderate explosion energies of about 1 x 
10^^ erg could be obtained for the reference model, i.e. 
10.8 Mq using the hybrid E0S2. The models with an 
early PNS collapse and hence the early onset of the ex- 
plosion reach smaller explosion energies. The shorter 
mass accretion period post bounce lead to less compact 
PNS configurations and therefore less steep density gra- 
dients at the PNS surfaces until they become gravitation- 
ally unstable and collapse due to the presence of quark 
matter in the interiors. The comparison of the two hy- 
brid EOSl and E0S2 {as = 0) is shown in Fig. [JT] for 
the 10.8 Mq progenitor model. The PNS obtained for 
B^/-* = 165 MeV (E0S2) is more compact, indicated 
by the more massive quark core in Fig. [17] (f ) of about 
0.15 Mq at the moment of shock breakout. It results in 
the higher density and temperature shown in the Figs.lTTl 
(e) and (b). Furthermore, the longer mass accretion time 
for the model using B^/^ = 165 MeV (E0S2) relates to 
lower densities and a steeper density gradient, surround- 
ing the central quark core of the PNS (see Fig.[T7](b)) at 
the moment of shock breakout. This allows for a stronger 
shock acceleration of the second shock wave expanding 
along the decreasing density gradient at the PNS sur- 
face. These differences are the origin of the higher ex- 
plosion energy estimates using E0S2 in comparison to 
EOSl. Hence, the acceleration of the formed hydrody- 
namic shock at the PNS surface is less intense for the 
models using the low bag constant. Lower explosion en- 
ergies are also obtained for the more massive progenitors. 

* Note that the values of pc, Tc and Ye are not the critical 
conditions for the appearance of quark matter. 
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Fig. 16. — Selected radial profiles for the 10.8 Mq reference model 
Illustrating is the back and forth propagating accretion shock at the 

The more massive envelopes of the PNSs with signifi- 
cantly higher densities prevents a similar strong shock 
acceleration. In general, the explosion energy estimates 
may be lifted to larger values if multi-dimensional phe- 
nomena (e.g. convection and rotation) are taken into 
account. 

4.2. PNS collapse to a black hole 

A special model is the 15 Mq progenitor, using the hy- 
brid E0S2 (see Table [2|). The evolution of this model is 
shown in Fig. [18] via radial profiles of velocity, density and 
temperature in the graphs (a), (b) and (c) respectively 
at four selected post bounce times. The graphs illustrate 
the last stable configuration as well as the PNS collapse. 
The evolutionary scenario, i.e. the PNS collapse due to 
the softening of the EoS in the mixed phase and the sub- 
sequent shock formation due to the stiffening of the EoS 
in the pure quark phase (clearly identified via the veloc- 
ity, density and temperature profiles in Fig. [18] (a), (b) 
and (c) respectively at the first three post bounce times), 
is in qualitative agreement with the reference case dis- 
cussed in § 3. During the shock propagation towards the 
PNS surface, mass accretion onto the second shock wave 
continuously increases the central density and hence the 
quark core of the PNS grows in mass. When the maxi- 
mum stable mass (given by the hybrid EoS) of the con- 
figuration is reached, the PNS collapses a second time as 
shown in Fig. [18] (dotted- lines). Relativistic effects be- 
come more and more important. This is illustrated in 




10 10 
Radius [i<m 

during the ongoing explosion phase at three late times post bounce. 
PNS surface established via to neutrino cooling/heating. 

Fig. [TS] (d) via the metric function a(t,a) = gu, i-e. the 
lapse function, which approaches zero. It indicates the 
appearance of the event horizon and hence the formation 
of a black hole. By our choice of a co-moving coordinate 
system, stable solutions for the evolution equations of 
energy and momentum cannot be obtained. Hence, the 
simulations cannot be continued beyond that point. 

The evolution scenario of black hole formation induced 
via the PNS collapse found via the hybrid EoS here is 
in qualitative agreeme nt wi th the results obtained by 
INakazato et"!!] (|2008b[) and INakazato et"!!] (|2010bf ) for 
very ma ssive stars of 100 Mp and for the 40 M© progeni- 
tor from lWooslev k, Weaveii ( )1995[ ). Furthermore, the re- 
sults agree qualitatively with the pure hadronic scenario 
of bla ck h ole formation discussed by ISumivoshi et al.l 
(|2007[) and lFischer et al.l (|2009D for several massive pro- 
genitor stars in the mass range of 40-50 M© . Differences 
occur in the larger maximum masses for the hadronic 
EoS. 

4.3. Robustness of the observable features 

The peak-like rise of the luminosities and mean neu- 
trino energies is an indirect consequence of the quark- 
hadron phase transition as discussed for the reference 
case in § 3. It was found to occur for all models under 
investigation. The results are shown in Figs. [T^] and [20] 
for the models using the hybrid EOSl (in graphs (a)) and 
E0S2 (in graph (b)). The different post bounce times 
(as listed in Table [3] together with the corresponding cen- 
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Fig. 17. — Selected radial profiles for the 10.8 Mq progenitor 
model, comparing the hybrid EOSl (red lines, at 248.57 ms post 
bounce) and EOS2 (blue linos, at 429.76 ms post bounce), both 
at the moments just before shock breakout along the decreasing 
density at the PNS surfaces. Graph (f) shows the composition, 
i.e. quark volume fraction (solid lines) and the hadronic matter 
(dashed lines). In addition to the radial profile in graph (f), the 
inset shows the composition with respect to the baryon mass. 

tral maximum densities Pmax) for the appearance of the 
second peak and the magnitude of the peaks depend on 
the microphysics and the progenitor modeL The quantity 
Pmax relates to the maximum central density obtained at 
the end of the PNS collapse. 

For all models under investigation (see Fig.lT^and l^O)) . 
the luminosities and mean neutrino energies follow a sim- 
ilar behavior. After a sharp rise on timescales on the or- 
der of milliseconds, the luminosities and mean neutrino 
energies increase and decrease in accordance to the ex- 
panding and descending standing accretion shock at the 
PNS surface, between 10-50 ms after the burst. A quasi- 
stationary state is obtained after about 100 ms after the 
sharp rise. The luminosities and mean neutrino energies 
decrease continuously indicating the ongoing explosions. 
All models evolve in a similar fashion as discussed in 
§ 3, through the PNS collapse and the subsequent evo- 
lution. Hence, the second peaks are dominated by 
and (t'^/r ; i^i//t) for ^ models under investigation (see 
Fig. I19p . Furthermore, the bursts are accompanied by 
significant increase of the mean neutrino energies (see 
Fig. [201). 

The magnitude and the width of the luminosities and 



mean neutrino energies of the second neutrino burst 
found are on the same order for all models under in- 
vestigation, independent of the progenitor model and in- 
dependent of the bag constant. On the other hand, the 
onsets for the second neutrino bursts take place at differ- 
ent times post bounce and depend on both the progenitor 
model and the bag constant. As shown in the Figs. [T51 
and[20l the second neutrino bursts for the 10.8 Mq and 
13 M0 models using EOSl occur at similar times because 
the conditions for the PNSs to become gravitationally 
unstable are obtained at similar post bounce times (see 
Table [2]). Since the structure of these two progenitor 
models is relatively similar, they evolve in a similar fash- 
ion and on a similar timescale post bounce until the PNSs 
collapse. The 15 M0 progenitor model on the other hand 
is more compact and the conditions for the PNS collapse 
using EOSl are obtained slightly earlier by about 40 ms 
(see Table [2] and compare with the post bounce time of 
the second neutrino burst in Figs. [TOlandEOl (a)). 

The simulations using the hybrid E0S2 is shown in 
the Figs, m (b) and [20] (b). They bring the differences 
between the progenitor models to light, due to the dif- 
ferent post bounce accretion times before the PNSs be- 
come gravitationally unstable and collapse. The lower 
progenitor masses relate to longer accretion times and 
hence later second neutrino bursts. The delay for the 
different models explored in this article is on the order 
of about 100 ms for the 10.8 to the 13 as well as between 
the 13 and 15 Mq progenitor models. In other words, 
the differences between the post bounce times for the re- 
lease of the second neutrino bursts relate to the same 
critical conditions for the quark-hadron phase transition 
but different evolutionary scenarios due to the different 
progenitor models. The future observation of such multi- 
peaked neutrino spectra might allow us to extract funda- 
mental information about the state of matter at extreme 
conditions. The magnitude of the second burst and its 
delay after the deleptonization burst after core bounce 
contains information about the quark and hadron EoSs. 
If the progenitor model is known and the pure hadronic 
EoS is fixed, it may be possible to extract the critical 
conditions for the quark-hadron phase transition from 
the observed neutrino signal. The post bounce times 
for the release of the second neutrino burst and the cor- 
responding central conditions for the PNS collapse are 
listed in Table [2] for all models under investigation. First 
results have been analyzed with respect to the possible 
observ ation of such a neutrino burs t bv iDaseupta et al.l 
()2010[) based on results discussed in lSagert et al.l ()2009D . 
exploring the 10.8 Mq progenitor model where the hybrid 
EOSl was used. They find that the operating neutrino 
detectors, Super-Kamiokande and IceCube, could detect 
such a z/g-burst from a future Galactic event. 

4.4. Corrections of the strong coupling constant 

Massive stars of 15 M© are in the expected mass range 
of exploding progenitors, that result in stable PNSs at 
least during the time of explosion. The evolutionary sce- 
nario of black hole formation applying the hybrid E0S2 
as discussed above in § 4.2, indicates that important in- 
gredients are missing. For that particular model, the 
PNS collapses to a black hole because the maximum sta- 
ble mass of the configuration, which is given by the hy- 
brid EoS, is reached during the collapse. The collapse in 
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Fig. 18. — PNS collapse of the 15 Mq progenitor model using EOS2, at selected post bounce times during the PNS collapse, shock 
formation and shock propagation as well as the black hole formation. 



turn is triggered due to the presence of quark matter in 
the interior. The inclusion of a^-corrections in the hy- 
brid EoS as introduced in §.2.3 allows for more massive 
hybrid star configurations. It increases the maximum 
stable mass significantly. Applying the hybrid E0S3 
— 155 McV, as — 0.3), leads to similar critical 
conditions for the onset of quark matter in comparison 
to E0S2, however with slightly higher critical density 
(see Table [5]). Furthermore, the post bounce time for 
the PNS to become gravitationally unstable and collapse 
(see Table [2]) is slightly delayed by about 50 ms (see 
Fig. [2T|) . which results in slightly higher temperatures 
(see Table [2]). The consequent evolutionary scenario is 
in qualitative agreement with the reference model as dis- 
cussed in §. 3, i.e. the 10.8 Mq using E0S2. The remain- 
ing hybrid star has a baryon mass of about 1.65 Mq at 
the moment when the simulation is stopped. Due to the 



larger maximum stable mass of the configuration, the 
PNS is found to be stable and does not collapse to a 
black hole. 

The evolution of the corresponding neutrino luminosi- 
ties and mean energies arc shown in Fig. 1211 which is 
in qualitative agreement with the spectra discussed for 
example of the reference model in §3. 

4.5. PNS structure and the early evolution 

The rise of the central density at the final phase of the 
Fe-core collapse of massive progenitor stars is shown in 
Fig. [22] (a) for the 10.8 Mq reference model. It reaches 
10^^ g/cm"^ at about 6 ms before bounce and increases 
rapidly on a timescale on the order of milliseconds up 
to nuclear densities, on the order of 2 — 4 x 10"'^^ g/cm^. 
In general, the maximum values obtained for the cen- 
tral density depend on the EoS and on the progenitor 
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Fig. 19. — Neutrino luminosities for the different progenitor models under investigation. The time gauges tc correspond to the moment 
when the maximum central densities are obtained for the different models during the PNS collapses. 



TABLE 3 

Estimated neutron star masses Mns and explosion energies 

Eexpl and the maximum CENTRAL DENSITIES pi,Ax FOR ALL 
MODELS UNDER INVESTIGATION. 



Prog. 


EoS 




^cxpl 




pmax 


[M©] 




[Me] 


[10^1 erg] 


[ms] 


[1015 












g/cm^] 


10.8 


EOSl 


1.431 


0.373 


248.78 


1.291 


10.8 


EOS2 


1.479 


1.194 


429.81 


1.806 


13 


EOSl 


1.465 


0.232 


241.20 


1.323 


13 


EOS2 


1.496 


0.635 


364.06 


1.788 


15 


EOSl 


1.608 


0.420 


175.07 


1.487 


15 


EOS2 


1.641 


unknown 


277.10^ 


15.362 _^ 


15 


EOS3 


1.674 


0.458 ~ 


312.99" 


1.342~ 



''Neutron star mass (baryon mass) and explosion energy, esti- 
mated at several 100 ms after the onset of explosion. 

'^Post bounce times tc, correspond to the moments when the max- 
imum central densities pmax are reached at the end of collapse 
black hole formation before positive explosion energy is achieved 

"^time of black hole formation 

"^central density at the time of black hole formation 

model. The moment of bounce, i.e. t = 0, is deter- 
mined when the maximum central density is obtained 
(see Fig. [22] (a)). After bounce, the central density de- 
creases slightly on the same short timescale on the order 
of milliseconds, until after about 10 ms post bounce the 
central density increases back again. It relates to the 
stalling of the bounce shock. After bounce, the central 
density stays above nuclear saturation density (see the 



horizontal dotted line in Fig. [22]) . During the later post 
bounce evolution, which is determined by the presence of 
the standing accretion shock and hence mass accretion, 
the central density increases up to 6.3 x 10"'^^ g/cvc? at 
about 428 ms post bounce (see Fig. [22] (b)). At about 
428.5357 ms post bounce, the PNS contraction proceeds 
into a collapse. 

The density increases rapidly up to 1.29 x 10^^ g/cm"^ 
(see Fig. [21] (b)). The post bounce time when the max- 
imum central density is reached during the PNS col- 
lapse is labeled tc = 429.3 ms. After ic, the central 
density decreases again slightly, which relates to the ini- 
tial expansion of the second shock wave formed inside 
the PNS. At about 0.1 ms after tc, the central density 
starts to oscillate at a short timescale on the order of 
milliseconds, shown in Fig. [H] (b). This oscillation is 
due to the oscillating PNS interior, which is triggered by 
the mass accretion onto the outwards propagating accre- 
tion shock where nucleons are dissociated into quarks. 
Hence, the quark core continues to grow in mass. About 
5 ms after tc, the saturation value is obtained, where 
Pccntrai — 1-25 X 10"'^^ g/cm'^ which remains constant with 
respect to time on the order of 100 ms. The central den- 
sity will continue to rise only at later times and at a 
longer timescale, which is given by fall-back of material 
enclosed inside the mass cut and due to deleptonization 
(i.e. coohng). 

The scenario explored here might reveal a connection 
between magnetars and neutron stars with quark matter 
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Mean neutrino energies for the different progenitor models under investigation. The same eonfiguration as Fig. 1191 



Fig. 

in the interiors. Due to flux conservation, the foUowing 
relation between the baryon density p and the magnetic 
field strength B applies, 



This expression indicates a B-field increase during the 
Fe-core collapse post bounce as well as an additional in- 
field increase during the quark-hadron phase transition, 
in accordance with the density increase. The evolution 
of the central magnetic field, according to the above re- 
lation, is shown in Fig. [53] with respect to time after 
bounce. It corresponds to the density evolution shown 
in Fig. [22l Illustrated are three different initial mag- 
netic field strengths, i.e. Bq = 10^ G, Ba = 10^° G and 
Bq = 10^^ G. These are values which m ight be obtained 
at pre-collapse stehar m odels (see e.g. IHeger fc Langeii 
l2000t iHegereFall [20051) . At core bounce, for these ini- 
tial configurations and according to the density evolu- 
tion as discussed above, the central B-fields reach values 
of 10^^, 10^"^ and 10^'^ G. The post bounce compression 
increases the B-fields additionally by a factor of about 
2 until the PNS collapses at about 428.3757 ms post 
bounce where, due to the density jump, the i3-fields rise 
to about 2.5 x 10^^"^'* G on a short timescale on the order 
of milliseconds. After the quark-hadron phase transition, 
the central density increases on a longer timescale on the 



order of 100 ms up to seconds and hence the magnetic 
field strength increases slowly on the same timescale. 

In addition to the evolution of the PNS interior in 
Fig.[23]for the central density evolution, close to the PNS 
surface the density decreases over several orders of mag- 
nitude. Hence we expect the magnetic field strengths at 
the PNS surface to be smaller in comparison to the evolu- 
tion at the PNS center. Surprisingly, the i?-fields reach 
values only slightly below 10^^~^^ G close to the PNS 
surface. During the later PNS evolution after the explo- 
sion has been launched, the PNS continues to contract. 
The subsequent density increase due to deleptonization 
on a longer timescale on the order of seconds, may cause 
an additional increase of the _B-field where based on the 
three initial choices used values of 10^^ G might be ob- 
tained at approximate simulation times on the order of 
10 seconds post bounce. 

In this ad-hoc approach, important effects such as the 
presence of rotation and fiuid instabilities have been ne- 
glected. Due to angular momentum conservation, the 
PNS spins up during the (second) collapse. Even if at 
the onset of the (second) collapse the PNS is in solid- 
body rotation, the fact that the collapse proceeds non- 
homologously leads to differential rotation. Since a ro- 
tating body has the least rotational energy when it ro- 
tates solidly, there is some free energy stored in dif- 
ferential rotation. The free energy in differential rota- 
tion may be converted into magnetic energy by winding 



Supernova explosions triggered by a quark-hadron phase transition 



25 



5r 

4.5 

_ 4 

0) 

CO o 
O 

^2.5 

I 1.5 
1 

0.5 

2- 



,1.5- 



i 0.5 




50- 




Fig. 21. — Luminosities and mean neutrino energies for the 
15 Mq progenitor model using the hybrid EOS3. The time gauge 
tc is again the moment when the maximum eentral density is ob- 
tained during the PNS collapse. 



up poloidal field into toroidal field. Further, depend- 
ing on the angular velocity profile, the PNS may also 
be s ubject to the magneto- rotational instability (MRI) 
(see iBalbus fc Hawlevlll998D . While field winding leads 
to a linear field strength growth, the MRI may lead 
to exponential field growth. Also convective motion 
from negative entropy and/or lepton gradients gener- 
ated by the propagation of the second shock and the sec- 
ond neutrino burst may lea d to dynamo action (see e.g. 
[Thompson fc DuncanI Il993[) . Therefore our estimated 
values for the evolution of the magnetic field strength 
represent a lower limit. In order to obtain a more pre- 
cise prediction of the magnetic field strength, detailed 
multidimensional investigations are necessary. 

Finally, Fig. [M] schematically shows the resulting PNS 
structure, with quark matter in the interior. More pre- 
cisely, the PNS is composed of an extended quark core 
surrounded by a high-density (including a reasonable 
amount of matter close to nuclear saturation density) 
hadronic mantle as well as a low-density envelope. The 
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Fig. 22. — Evolution of the central density for the 10.8 M@ ref- 
erence model with respect to time after bounce. The dotted line 
marks nuclear saturation density. 
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Fig. 23. — Evolution of the central magnetic field strength for 
the 10.8 Mq reference model with respect to time after bounce 
for three different initial magnetic field strengths (10^ G, 10^" G, 
10" G). 

surface of the PNS can be defined via the contracting 
ncutrinospheres. The evolution of the PNS is determined 
by fall-back of very low-density hadronic material en- 
closed inside the mass cut. The mass cut estimate is a 
dynamical quantity which changes during the evolution. 
The mass enclosed inside the mass cut will form the rem- 
nant hybrid star. The remaining hybrid star cools after 
mass accretion vanishes, initially via deleptonization up 
to about one minute and later via photons up to millions 
of years. Matter outside the mass cut has become grav- 
itationally unbound during the ongoing explosion and 
mass ejection from the PNS surface via neutrino heating 
(i.e. the neutrino-driven wind). This mass will be ejected 
into the interstellar medium. As can be seen from Fig. l24[ 
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Fig. 24. — Schematic composition, illustrating the radial profile 
of the remaining hybrid PNS configuration for the reference 10.8 
M0 progenitor model at late post bounce time. 

a reasonable amount of mass (about 7 x 10~^ M0) which 
corresponds to the ejected matter, has become neutron- 
rich with an electron fraction of Yf, ~ 0.3. During the 
later expansion, the electron fraction increases and val- 
ues between ~ 0.35 — 0.5 are obtained. Detailed explo- 
sive nucleosynthesis investigations are required in order 
to confirm or exclude the scenario as a possible site for 
the synthesis of heavy elements via the r-process. 

4.6. The explosion ejecta 

The mass elements which are ejected directly (ballisti- 
cally). are located at several 10^ km during the Fe-corc 
collapse, bounce and until the quark-hadron phase tran- 
sition. They have densities and temperatures of several 
10^ g/cm^ and about 0.2 MeV, which remain constant 
with respect to time. These mass elements are relatively 
far away from the center and belong to the surrounding 
part of the Fe-core. The composition is dominated by 
^®Si and ^^S, where Ye ~ 0.5. The structure of these 
layers is given by the progenitor model and changes only 
insignificantly during the collapse, bounce and early post 
bounce phases. 

The situation changes when the expanding explosion 
shock crosses these mass elements. This happens at a 
post bounce time of about 450 ms for the 10.8 Mq ref- 
erence model. The mass elements contract rapidly to- 
wards the PNS surface to about 10 km, during the PNS 
collapse. Density and temperature increase to several 
10^^ g/cnT' and several tens of MeV, on a short timescale 
on the order of milliseconds. Furthermore, the entropy 
per baryon increases from values below 10 k^ to about 
60 ks and the electron fraction decreases from Yg ~ 0.5 
to Ye ~ 0.1. During the later expansion, density and tem- 
perature decrease initially on the same short timescale. 



where a larger value for the electron fraction is obtained, 
between Ye ~ 0.33 — 0.48, and the entropy reduces to val- 
ues close to 40 ks per baryon. After the initial rapid ex- 
pansion, these direct explosion ejecta expand on a longer 
timescale on the order of 100 ms up to seconds. Den- 
sity and temperature decrease on the same timescale and 
reach values between 10'^"''' g/cm"^ and below 0.1 MeV at 
about 2 seconds post bounce. For these obtained condi- 
tions, the reaction rates become much smaller and hence 
the electron fraction changes on a longer timescale on the 
order of seconds. In addition, since the later expansion 
proceeds adiabatically, the entropy per baryon remains 
constant with respect to time as well, where values of 
~ 40 ke per baryon are obtained. 

Mass elements which are not ejected directly in the 
initial shock expansion fall back onto the PNS surface 
at the center. During the short infall timescale on 
the order of milliseconds, density and temperature in- 
crease up to 10^"^ g/cm*^ and 5-10 MeV. During the 
mass accretion onto the PNS surface, weak-equilibrium 
establishes at a low value of the electron fraction of 
Ye ~ 0.05 and the entropy per baryon decreases to val- 
ues below 5 kg. However, after about 1.5 seconds post 
bounce t he neutrino-dr i ven w in d appears as d iscusse d 
in e.g. lArcones et all (I2007D : iFischer et all (|2010bD . 
One of the most important result of the appearance 
of the neutrino-driven wind relates to the high entropy 
conditions obtained for the accelerated matter. It re- 
sults in gene rally proton-rich m a terial , where Ye ~ 0.55 
(first seen inlLiebendorfer et al.l (12003 ) and e xplained in 
iFrohhch et all (|2006[) and lPruet at all (|2006[ )). 



5. SUMMARY 

We performed core-collapse supernova simulations of 
massive stars in the mass range of 10.8 to 15 Mq, 
based on general relativistic radiation hydrodynamics 
employing a sophisticated EoS and three flavor Boltz- 
mann neutrino transport. A description for three flavor 
quark matter is implemented based on the bag model 
with the different parameter choices, {B^^^ = 162 MeV, 
a, = MeV), {B^/* = 165 MeV, a, = MeV) and 
(^1/4 = 155 MeV, as = 3 MeV). The resulting quark 
EoSs, EO Sl, E0S2 and E0 S3, are coupled to the hadron 
EoS from lShen et al.l (|1998f ). where Gibbs conditions are 



applied for the transition between hadron and quark mat- 
ter. It results in an extended co-existence region in the 
phase diagram where nucleons and quarks are present, 
the mixed phase. The thermodynamic conditions for the 
appearance of quark matter relate to critical densities 
close to nuclear matter density, for temperatures of tens 
of McV and a low proton-to-baryon ratio. 

During the evolution of massive stars, we find that 
quark matter appears already at the Fe-core bounce for 
the chosen parameters. However, the quark matter vol- 
ume fraction is found to be rather small with x ^ 0.1. 
The quark matter volume fraction increases slowly over 
several 100 ms during the early post bounce mass accre- 
tion phase, during which central density and temperature 
increase while the central electron fraction decreases. 
The EoS in the mixed phase is significantly softer com- 
pared to the pure hadronic and the pure quark phases. 
Hence, the PNS contraction accelerates during the post 
bounce evolution as more and more matter of the PNS in- 
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terior is converted into the mixed phase. Finally, the cen- 
tral PNS configuration becomes gravitationally unstable 
and the contraction proceeds into an adiabatic collapse. 
Density and temperature increase on a short timescale 
on the order of milliseconds, which in turn favors quark 
matter over hadronic matter. The stiffening of the EoS 
in the pure quark phase halts the collapse and a strong 
hydrodynamic shock wave forms. The shock appears ini- 
tially as a pure accretion front. It propagates outward 
along the decreasing density of the PNS. The propaga- 
tion is determined via the balance of ram pressure from 
the supersonically infalling material ahead of the shock 
and the thermal pressure of the quarks as well as neutrino 
heating behind the shock. At the PNS surface where the 
density decreases over several orders of magnitude, the 
accretion shock accelerates and positive matter velocities 
are obtained. The accretion front turns into a dynamic 
shock wave, which continues to expand. This moment 
establishes the onset of explosion. It finally merges with 
the standing accretion shock from the Fe-core bounce, 
which remained unaffected from the happenings inside 
the PNS. All models under investigation follow the same 
evolutionary behavior and lead to explosions, if the PNS 
mass does not exceed the maximum stable mass of the 
configuration given by the hybrid EoS. 

Investigating the large uncertainty in core-collapse su- 
pernova input physics, i.e. the state of matter at high 
densities and temperatures and with a low proton-to- 
baryon ratio, we were able to explore a new explosion 
mechanism in simulations of massive stars. It serves 
as an addition to the well explored neutrino-driven, the 
magnetically-driven and the acoustic mechanisms. Note 
that these appear to be working only in multiple spatial 
dimensions (except for the low mass 0-Ne-Mg-core) , due 
to the nature of the physics directly involved in these ex- 
plosion mechanisms. The explosions obtained by taking 
QCD degrees of freedom into account are presumably not 
restricted to dimensional limitations. However, moderate 
explosion energies on the order of 10^^ erg could only be 
obtained for the 10.8 M0 progenitor model and for a par- 
ticular choice of parameters. Taking multi-dimensional 
phenomena into account, such as rotation and the de- 
velopment of fiuid instabilities, may shift the explosion 
energies to robust values. 

The sudden density rise during the seconds collapse 
indicates an additional increase of the magnetic field 
strength. It enables us to give a lower estimate of the 
expected magnetic field strength for protoneutron stars 
with quark matter cores. However, with initial magnetic 
fields that are confirmed by stellar evolution calculations, 
a connection between hybrid stars and the special neu- 
tron star class Magnetars with magnetic fields up to sev- 
eral 10"'^^ G could not be established. Important aspects 
of the magneto-hydrodynamics evolution have not been 
taken into account due to the current restriction to spher- 
ical symmetry. They will be explored in a detailed study 
in an upcoming article. 

Direct observables which allow us to identify the ap- 
pearance of quark matter in the PNS interior cannot be 
expected in the neutrino signal, because matter is opaque 
for neutrinos. Neutrinos can only diffuse out of the PNS 
interior on timescales on the order of seconds. How- 
ever, the presence of a strong hydrodynamic shock wave 
changes the situation. It releases an additional outburst 



of neutrinos when crossing the neutrinospheres. This 
millisecond neutrino burst is dominated by i^e, because 
matter is neutron-rich where positron captures dominate 
over electron captures, and by similar fluxes of Ve and 
[v^/t,v^/t) that are produced via pair process. The de- 
lay of the second neutrino burst to the deleptonization 
burst from the Fe-core bounce and its magnitude, contain 
information about the progenitor model, the hadronic 
and quark EoSs, including the critical conditions for the 
onset of quark matter as well as the nature of the phase 
transition from hadronic matter to quark matter. It will 
become observabl e for future Ga l actic events, as con- 
firmed recently bv lDasgupta et all ()2010D for the operat- 
ing neutrino detectors Super-Kamiokande and IceGube. 
The knowledge of the progenitor model and the hadronic 
EoS may allow us to decode hidden information about 
the state of matter at high densities and temperatures, 
conditions which will likely be reached in future heavy- 
ion collision experiments at FAIR/GSI (Germany) and 
NIGA/Dubna (Russia). Furthermore, direct observables 
from the quark-hadron phase transition can be expected 
in the emission of gravitational waves. These might be- 
come observable in the near future if the sensitivity im- 
provements of the gravitational wave detectors continue. 
They will serve as important addition to the neutrino 
spectra. 

The composition of the ejecta can also be used to probe 
the explosion scenario. Although core-collapse super- 
novae have been considered as the favored site for pos- 
sible r-process nucleosynthesis for a long time, the re- 
quired conditions have not been found in simulations so 
far. Therefore it would be interesting to explore whether 
an r-process can ensue in the layers being ejected here. 
We defer, however, a thorough investigation to a sep- 
arate paper. The reduced nuclear reaction network in- 
cluded here to account for energy generation is not suited 
to follow the Ye evolution and nucleosynthesis at late 
times. Nevertheless, there is only a small range of lay- 
ers to be considered. Neutron-rich matter from the deep 
layers close to the PNS surface attains ~0. 33-0. 48 af- 
ter shock heating and is ejected ballistically. Layers in 
a second class are not directly ejected, rather they fall 
back onto the PNS surface, where they experience con- 
tinued neutrino heating and can subsequently be ejected. 
Initially, these layers can be very neutron-rich but turn 
proton-rich as they expand in the neutrino-wind devel- 
oping at later times. Therefore, only the directly ejected 
zones may contain viable conditions for neutron-rich nu- 
cleosynthesis, provided they are located not too far above 
the PNS surface. 

With this study we address a deeper understanding of 
the inclusion of QGD degrees of freedom in radiation hy- 
drodynamics studies of astrophysical scenarios, explored 
at the example of core-collapse supernova simulations of 
massive stars. The quark matter description applied here 
leads to an early onset of the chiral phase transition. 
This is an active subject of research where several as- 
pects, such as the chiral phase transition in general, the 
hadronic fre eze out, the existence of a possible quarky- 
onic phase (jMcLerran &: Pisarskil [2007[ ) and the nature 
of the transition between hadronic and quark matter, 
are investigated. Further developments of quark matter 
EoSs where different conditions with respect to the crit- 
ical conditions for the onset of deconfinement may favor 
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different evolutionary scenarios, will be studied in the 
future. Furthermore, simulations in multiple spatial di- 
mensions are required in order to support the findings of 
the present article. 
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